ODC  FILE  COE#  AD  AO  65831 


S FWF 


3Bk.ja 


D D C 
j7"> 'xSIrT 

i mar  15  1979  i 


■:  MAR 


1 his  dor.unv.-a? 
for  public  • i 

distribution  is  ui 


:n  ( 


_ T THERMAL-MECHANICAL  AND  THERMAL 

' - — i ■»_  — 

J5EHAVIOR  OF  HIGH-TEMPERATURE  STRUCTURAL  ^MATERIALS 


Interim  j/eplwrt,  fcorCWR  ^ 


pec.  9^ 


Agreement  NOj^-J^NppQlA^S-C-fiASl) ' 

Q$)- 

April  1,  197o^-  December  31,  1978 


&T7.  p.„, 


y Hasselman  and  M.P./Kamat  j jp|  *7 


— it»j«oop«-ration  with 

E.K.j Beauchamp,  L./Bentsen^K./chyung|  J.  Heinrich,  J.W.  McCauley, 
-ffTTCnoch,  KT- SaTyaMUrt liyT'JTPT'sI n gTi , J.C.  Swearengen,  J.R.  Thomas , 
G.E.  Youngblood,  Y.  Tree  and  W.A.  Zdaniewski  . 


Departments  of  Materials  Engineering  and 
/•'“  Engineering -'ssience  and  Mechanics 

7 Virginia  Polytechnic  Institute  and  State  University^) 
^ — __  Blacksburg,  Virginia  24061 


-4U1  7 


■ pO  7 70t> 


\ p 


SECIilUl  < CL  »SSI  FIC  ATIOH  OF  THIS  PAGE  (HTi«n  Dala  Entered) 


REPORT  DOCUMENTATION  PAGE 


. REPORT  N 'MISER 


4.  TITLE  5uftfllf<0 


12.  GOVT  ACCESSION  NO. 


RF.AD  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


. RECIPIENT'S  CATALOG  NUMflER 


Thermo-Mechanical  and  Thermal  Behavior  of  High- 
Temperature  Structural  Materials 


7.  AUTHOR!*} 

D.  P.  H.  Hasselman  and  M.  P.  Kamat 


9.  performing  organization  name  and  address 


6 PERFORMING  ORG  REPORT  NUMBER 


9.  CONTRACT  OR  GRANT  NUMSERfl) 


10  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A WORK  UNIT  NUMBERS 


Virginia  Polytechnic  Institute  and  State  Universit) 
Blacksburg,  VA  24061 


M.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


12.  REPORT  DATE 

December  31,  1978 


13.  NUMBER  OP  PAGES 


I 


MONITORING  \GFNCY  NAME  A AbtfRESS(irt/lfT*rcnl  from  Controlling  Olt Ico)  15.  SECURITY  CLASS,  (ol  Ihlo  roport) 

Unclassified 


IS*.  DECL  ASSIFICATION' DOWNGRADING 
SCHEDULE 


H DISTRIBUTION  STATEMENT  (ol  I M»  Report) 


See  distribution  list. 


17  DISTRIBUTION  STATEMENT  (ol  the  ebelrecl  entered  In  Block  30,  II  dlllerenl  Ir on  Report) 


19.  KEY  WOROS  (Continue  on  rereree  aid*  II  nacaaaary  and  Identity  by  block  nuaiSorJ 

Absorption  coefficient,  micro-cracking,  thermal  stress  resistance,  thermal  fati- 
gue, figure-of-merit , ice,  thermal  conductivity,  thermal  diffucivity,  glass- 
ceramic,  mica-alumina,  borosilicate  glass,  thermal  buckling,  reaction-sintered 
silicon  nitride 


20.  AOSTRACT  (Continue  on  r+vmrao  oldo  H n»c*»»mry  and  Identity  by  block  numbmr) 


-OVER- 


1473  EDITION  OF  I NOV  «5  IS  OBSOLETE 


Unclnusi  rii  d 

SECURII.Y-CLASSIFICATION  OF  THIS  PAGF  iMinn  no  la  Fnleredl 


dLUJ  HI  T Y CLASSIFICATION  OF  THIS  PAOEfl*h»n  Dmlt  Enfte  ■!) 


This  report  presents  the  results  obtained  over  the  period  of  4/1/78  to 
12/31/78  on  a research  program  on  the  thermo-mechanical  and  thermal  behavior 
of  high-temperature  structural  materials.  A total  of  eleven  studies  were  com- 
pleted and  reported  in  the  form  of  individual  chapters  as  follows: 

I.  "Thermal  Stress  Analysis  of  Partially  Absorbing  Brittle  Ceramics 
Subjected  to  Radiation  Heating". 

II.  "Analysis  of  Thermal  Stress  Resistance  of  Micro-Cracked  Brittle 
Ceramics". 

III.  "Role  of  Physical  Properties  in  the  Resistance  of  Brittle  Ceramics  to 
Fracture  in  Thermal  Buckling" 

IV.  "Thermal  Stress  Resistance  Parameters  for  Brittle  Materials  Subjected 
to  Thermal  Stress  Fatigue" 

V.  "On  the  Thermal  Fracture  of  Ice". 

VI.  "Effect  of  Cracks  on  Thermal  Conductivity", 

VII.  "Effect  of  Oxidation  on  Thermal  Diffusivity  of  Reaction-Sintered 
Silicon  Nitride". 

VIII.  "Effect  of  Crystallization  on  the  Thermal  Diffusivity  of  a Cordierite 
Glass-Ceramic", 

IX.  "Thermal  Diffusivity  of  Ba-Mica/Alumina  Composites". 

X.  "Effect  of  Alumina  Dispersions  on  the  Thermal  Conductivity/Dif fusi- 
vity  and  Thermal  Stress  Resistance  of  a Borosilicate  Glass", 

XI.  "Figures-of-Merit  for  the  Thermal  Stress  Resistance  of  High-Tempera- 
ture Brittle  Ceramics:  A Review" 


■ 


Unclassi fied 

SECURITY  CLASSIFICATION  of  THIS  RAOinWiFd  DfI* 


Ttn 


THERMO-MECHANICAL  AND  THERMAL 
BEHAVIOR  OF  HIGH-TEMPERATURE  STRUCTURAL 
MATERIALS 


PREFACE 

\ 

} j 

Technical  ceramics  because  of  their  chemical  inertness,  high  melting 
point,  good  wear  resistance,  excellent  mechanical  stability  at  high  tempera- 
ture and  other  unique  properties,  represent  a class  of  materials  eminently 
suited  for  many  critical  engineering  applications.  Unfortunately,  because 
of  their  brittleness  and  unfavorable  combination  of  pertinent  material  pro- 
perties, technical  ceramics  generally  are  highly  susceptible  to  catastrophic 
failure  in  non-uniform  thermal  environments,  which  give  rise  to  thermal 
stresses  of  high  magnitude. 

Thermal  stress  failure  analysis  of  structural  materials  represents  a 
multi-disciplinary  problem  which  involves  the  principles  of  heat  transfer , 
mechanics  and  materials  engineering.  Over  the  last  few  decades  much  general 
understanding  of  the  nature  of  thermal  stress  failure  of  brittle  materials 
has  been  generated.  However,  due  to  the  multi-disciplinary  nature  of  the 
problem,  the  ability  to  predict  thermal  stress  failure  quantitatively  for 
design  or  other  purposes  has  lagged  behind  the  progress  made  in  other  engi- 
neering fields.  The  objective  of  the  present  program  is  to  improve  the 
qualitative  and  quantitative  understanding  of  the  nature  of  thermal  stress 


failure  of  brittle  structural  materials.  In  order  to  achieve  this  objective, 
the  participating  investigators  and  scopC'''o¥ythe  program  are  organized  in  a 


manner  which  takes  full  advantage  of  the  combined  inputs  from  a number  of 
engineering  disciplines  including  mechanical  engineering,  applied  mechanics 
and  materials  science  and  engineering.  This  report  presents,  in  the  form 
of  individual  chapters,  completed  studies  prepared  for  publication  over  the 
contracting  period  of  4/1/78  to  12/31/78. 

The  effort  devoted  to  this  program  can  be  divided  into  three  separate 
individual  areas,  namely:  (A)  Analysis  of  thermal  stress  failure;  (B)  mea- 
surement and  analysis  of  material  properties  relevant  to  thermal  stress  failure 
and  (C)  reports  with  the  objective  of  the  dissemination  of  technical  informa- 
tion on  thermal  stress  failure.  The  completed  studies  are  grouped  together 
in  these  three  major  technical  areas  with  brief  descriptions  of  current  work 
in  progress  as  follows: 

A.  Analysis  of  Thermal  Stress  Failure. 

Chptr.  I.  D.  P.  H.  Hasselraan,  J.  R.  Thomas,  Jr.,  M.  P.  Kamat  and  K. 

Satyamurthy,  "Thermal  Stress  Analysis  of  Partially  Ab- 
sorbing Brittle  Ceramics  Subjected  to  Radiation  Heating," 

J.  Am.  Ceram.  Soc.,  (in  review). 

Chptr.  II.  D.P.H.  Hasselman  and  J.P.  Singh,  "Analysis  of  Thermal 
Stress  Resistance  of  Micro-cracked  Brittle  Ceramics," 

J.  Am.  Ceram.  Soc.,  (in  review). 

Chptr.  III.  D.P.H.  Hasselman,  "Role  of  Physical  Properties  in  the 
Resistance  of  Brittle  Ceramics  to  Fracture  in  Thermal 
Buckling,"  J.  Am.  Ceram.  Soc.,  (in  press). 

Note:  This  study  was  initiated  under  previous  support 
provided  by  ARO  and  revised  substantially  as  part  of  the 


present  program. 


Chptr.  IV.  D.P.H.  Hasselman  and  W.A.  Zdaniewski,  "Thermal  Stress 
Resistance  Parameters  for  Brittle  Materials  Subjected 
to  Thermal  Stress  Fatigue,"  J.  Am.  Ceram.  Soc.,  61  (7-8) 
375  (1978). 

Chptr.  V.  D.P.H.  Hasselman  and  Y.  Tree,  "On  the  Thermal  Fracture  of 
Ice,"  J.  Mat.  Sc.,  (in  press). 

Research  currently  underway  in  this  phase  of  the  program  includes  a 
number  of  studies  involving  finite  element  modelling  of  transient  thermal 
stresses  and  an  analysis  of  non-linear  convective  heat  transfer  in  thermal 
stress  failure  of  brittle  materials. 


B.  Measurement  and  Analysis  of  Material  Properties  Relevant  to  Thermal 
Stress  Failure. 


Chptr.  VI. 


Chptr.  VII. 


Chptr.  VIII. 


Chptr.  IX. 


Chptr.  X. 


D.P.H.  Hasselman,  "Effect  of  Cracks  on  Thermal  Conductivity," 

J.  Comp.  Mat.,  12,  403-07  (1978). 

W.  Zdaniewski,  H.  Knoch,  J.  Heinrich  and  D.P.H.  Hasselman, 
"Effect  of  Oxidation  on  Thermal  Diffusivity  of  Reaction- 
Sintered  Silicon  Nitride,"  Ceram.  Bull,  (in  press). 

K.  Chung,  G.E.  Youngblood  and  D.P.H.  Hasselman,  "Effect 

of  Crystallization  on  the  Thermal  Diffusivity  of  a Cordierite 
Glass-Ceramic,"  J.  Amer.  Ceram.  Soc.,  (in  press). 

G.E.  Youngblood,  L.  Bentsen,  J.W.  McCauley  and  D.P.H. 
Hasselman,  "Thermal  Diffusivity  of  Ba-Mica/Alumina  Com- 
posites," J.  Amer.  Ceram.  Soc.,  (in  press). 

D.P.H.  Hasselman,  J.C.  Swearengen,  E.K.  Beauchamp  and 
W.A.  Zdaniewski,  "Effect  of  Alumina  Dispersions  on  the 


Thermal  Conductivlty/Dif f usivity  and  Thermal  Stress  Re- 
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sistance  of  a Borosilicate  Glass,"  J.  Am.  Ceram.  Soc., 

(in  review). 

It  should  be  noted  that  the  above  studies  were  carried  out  in 
cooperation  with  other  organizations  and  also  received  considerable 
internal  VPI  support.  Chptrs.  VIII  and  IX  represent  completion  of 
brief  studies  carried  out  under  previous  ONR-support. 

As  part  of  the  on-going  research,  the  observations  reported  in 
Chptr.  VIII  are  being  pursued  further  in  the  form  of  quenching  studies 
of  reaction-sintered  silicon  nitride  subjected  to  various  oxidation 
treatments.  Experimental  data  obtained  at  MERDI  (under  previous  ONR- 
support)  on  a temperature  dependent  boundary  resistance  of  glass-metal 
composites  with  a rather  strong  positive  temperature  dependence  of  the 
thermal  diffusivity  are  being  analyzed.*  Also,  thermal  properties  of  a 
polycrystalline  alumina  used  in  experimental  thermal  stress  studies 
are  being  determined.  Other  measurements  in  progress  serve  to  verify 
the  analytical  findings  reported  in  Chptrs.  II  and  VI. 

C.  Dissemination  of  Technical  Information  on  Thermal  Stress  Fracture  of 
Brittle  Materials. 

Chptr.  XI.  D.P.H.  Hasselman,  "Figures-of-Merit  for  the  Thermal 

Stress  Resistance  of  High-Temperature  Brittle  Materials," 
Ceramurgia  International  (in  press) 

* Of  interest  to  note  is  that  similar  observations  are  being  made  in  parallel 
studies  on  potential  MHD-materials : J.J.  Rasmussen,  Task  N,  "Materials  Evalu- 
ation," Report  to  DOE-MHD,  Jan. -March,  1978. 


*L-  £ 


Effort  currently  underway  as  part  of  this  phase  of  the  program 
is  devoted  to  the  assembling  of  literature  references  and  other  informa- 
tion on  thermal  stress  fracture  of  brittle  materials  for  the  purpose  of 
compiling  an  extensive  bibliography  and  the  preparation  of  an  in-depth 
review  acticle.  Furthermore,  a conference  is  being  organized  with  the 
tentative  title:  "Thermal  Stresses  in  Structures  and  Materials  in  High 
Heat  Flux  Environments."  This  conference  will  take  place  in  the  early 
part  of  1980.  Exact  dates  and  details  of  the  program  will  be  announced 
shortly. 
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CHAPTER  I. 

THERMAL  STRESS  ANALYSIS  OF  PARTIALLY  ABSORBING 
BRITTLE  CERAMICS  SUBJECTED  TO  RADIATION  HEATING 
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THERMAL  STRESS  ANALYSIS  OF  PARTIALLY  ABSORBING 
BRITTLE  CERAMICS  SUBJECTED  TO  RADIATION  HEATING 


by 


D.P.H.  Hasselman,  J.  R.  Thomas,  Jr.,  M.  P.  Kamat  and  K.  Satyamurthy 


Departments  of  Materials  Engineering,  Mechanical  Engineering 
and  Engineering  Science  and  Mechanics 
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ABSTRACT 


Thermal  stresses  in  partially  transparent  brittle  ceramics  sub- 
jected to  radiation  heat  transfer  were  analyzed  for  two  convective  cool- 
ing conditions.  The  absorption  coefficient  is  introduced  as  a material 
property  which  affects  thermal  stress  resistance.  Appropriate  thermal 
stress  resistance  parameters  were  derived. 


I.  INTRODUCTION 
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Brittle  ceramic  materials  for  structural  purposes  are  well-known  to 
be  susceptible  to  catastrophic  failure  under  conditions  involving 
thermal  stress^.  For  purposes  of  proper  design  and  selection  of  mater- 
ials for  structures  subjected  to  thermal  stress,  it  is  critical  that  the 
role  of  the  physical  properties  which  affect  thermal  stress  resistance 
of  brittle  materials  is  well  understood.  For  this  purpose,  literature 
solutions  were  presented  for  the  role  of  the  pertinent  physical  proper- 
ties which  control  thermal  stress  fracture  for  conditions  of  steady- 

2 3 4 

state  or  transient  convective  heat  transfer  , thermal  buckling  and 

• « 

crack  propagation^.  For  the  solutions  presented,  so-called  "thermal 
stress  resistance"  parameters  were  derived^’ ^ on  the  basis  of  which  the 
optimum  material  can  be  selected  for  prescribed  heat  transfer  conditions 
and  failure  criteria. 

At  sufficiently  high  temperatures  heat  transfer  by  radiation  also 
can  play  a significant  role  in  the  thermal  stress  fracture  of  brittle 
materials.  Previously,  analytical  solutions  were  presented  for  the 
thermal  stress  resistance  of  opaque  materials  for  which  the  incident 

g 

black-body  radiation  was  absorbed  in  the  immediate  surface  . A further 

solution  addressed  itself  to  a material  subjected  to  radiation  heating 

which  was  transparent  below  a prescribed  wave-length  and  opaque  above 

9 

this  value  of  wave-length  . As  far  as  the  authors  are  aware  to  date  no 
thermal  stress  analysis  has  been  presented  for  ceramic  materials  sub- 
jected to  radiant  energy  in  which  the  incident  radiation  is  being 
transmitted  and  absorbed  throughout  the  material.  Such  an  analytical 
solution  is  presented  in  this  paper,  with  primary  emphasis  on  establish- 
ing the  role  of  the  appropriate  physical  properties  which  determine  the 

1 

— J 


thermal  stress  fracture  behavior.  The  boundary  conditions  chosen  for 
the  analysis  are  those  which  will  meet  the  objectives  of  the  study  in 
the  most  simple  and  expedient  manner. 

II.  THEORY 

A.  Boundary  Conditions 

The  analysis  will  be  carried  out  for  an  infinite  flat  plate  of 
thickness  2a  located  in  the  yz-plane  and  -a  < x < a.  At  x = a,  -a  the 
plate  is  subjected  at  t = 0 to  a steady-state  perpendicularly  incident 
uniform  radiation  such  as  from  a laser  or  focused  black-body  radiation. 
The  optical  properties  of  the  material  will  be  considered  "grey",  i.e., 
independent  of  wave-length.  The  reflectivity  of  the  material  will  be 
assumed  sufficiently  low,  so  that  the  effect  of  multiple  internal  re- 
flections within  the  plate  on  the  absorbed  energy  can  be  neglected.  The 
choice  of  symmetric  heating  of  the  plate  simplifies  the  expressions  for 
the  thermal  stresses  as  the  moments  of  the  temperature  distributions 
around  x = 0 vanish  and  need  not  be  considered.  The  properties  which 
control  the  magnitude  of  thermal  stress,  such  as  the  coefficient  of 
thermal  expansion.  Young's  modulus,  Poisson's  ratio  and  the  thermal 
conductivity  and  diffusivity,  also  are  assumed  independent  of  tempera- 
ture. Other  boundary  conditions  or  wavelength  and/or  temperature-depen- 
dent materials  properties  can  be  considered  by  modifications  of  the 
present  analysis  or  by  numerical  methods. 

Throughout  the  analysis  it  will  be  assumed  that  the  temperature  of 
the  plate  is  sufficiently  low  so  that  the  radiant  energy  emitted  by  the 
material  is  very  small  compared  to  the  incident  radiation.  For  justifi- 
cation of  this  assumption,  it  should  be  noted  that  in  transient  heat 
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transfer,  generally  the  thermal  stresses  reach  their  maximum  values  dur- 
ing the  initial  stages  of  the  transient  temperature  history.  For  the 
derivation  of  thermal  stresses  then,  the  plate  can  be  considered  to  be 
subjected  to  a constant  heat  flux.  This  same  approach  was  taken  in  the 
analysis  of  thermal  stresses  in  opaque  bodies  heated  by  radiation  and 

g 

validated  by  numerical  examples  . Of  course,  for  an  estimate  of  the 
total  temperature  history  or  if  the  intensity  of  incident  radiation  is 
sufficiently  high  that  surface  melting  or  visco-elastic  effects  occur 
before  brittle  fracture  can  take  place  in  the  cooler  interior  of  the 
plate,  the  assumption  of  constant  radiant  heat  flux  for  thermal  stress 
analysis  will  require  modification. 


B.  Derivation  of  Transient  Temperatures. 

For  normally  incident  radiation  intensity  qQ  outside  of  the  plate, 
the  heat  absorbed  by  the  surface  equals  cqQ,  where  e is  the  emissivity 
or  E = (1  - r)  where  r is  the  reflectivity. 

The  intensity,  q,  of  the  radiation  after  entering  the  plate  at  x = 


a is: 

_ -y(a-x)  /1N 

q = £q  e (1) 

o 

where  y is  the  absorptivity. 

Similarly  the  intensity  due  to  radiation  entering  the  plate  at  x 
= -a  becomes: 

q = eqoe-y(a+x)  (2) 


The  rate  of  heat  absorption  g'"  at  x is  the  sum  of  the  heat  absorbed 
due  to  the  fluxes  of  Eqs.  (1)  and  (2)  and  is  expressed: 


in  „ -ya,  yx  . -yx. 
g = yeqQe  (e  + e H ) 


(3) 


which  can  be  written: 
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g"'  = 2peq  e ^acosh(yx) 


For  the  derivation  of  the  transient  temperatures,  the  following 
differential  equation  needs  to  be  solved^: 


3 T g"'  (x)  1 9T 

a 2 K 6 3t 

dX 


(5) 


where  T is  the  temperature,  K is  the  thermal  conductivity  and  6 is  the 
thermal  diffusivity. 

Solutions  were  obtained  for  two  limiting  cases  for  external  cooling 
conditions: 

1.  A constant  surface  temperature  at  x = a,  -a,  which  physically 
corresponds  to  the  plate  being  contained  in  a convective  heat 
transfer  environment  with  a heat  transfer  coefficient  h = 00 . 

2.  The  complete  absence  of  cooling  at  the  plate  surface  which 
corresponds  to  an  external  convective  heat  transfer  environment 
with  heat  transfer  coefficient  h = 0. 

For  these  two  cases,  the  transient  temperature  distributions  can  be 
derived  as  follows: 

1.  h = 00 

For  this  case  the  boundary  conditions  can  be  expressed: 

3T 


3x  (°’  " ° 

T (x,  0)  = Tq 

T(a,  0)  = T(-a,  0)  = T 


(6a) 

(6b) 

(6c) 


Substitution  of  Eq.  (4)  into  Eq.  (5)  and  solution  by  mathematical 
techniques  common  to  the  solution  of  transient  temperatures^  yields: 


4ucq  e H cosh(ya)  n=°°  , -BXniw  i\n  A 

T(x,t)  - To  +-^ — l I*1'8  V(2  ° — 1 (7) 

n=o  (y  +X  ) X 

n n 


** is Si 


-r- 


where  A = (n  + l/2)ir/a 
n 


(8) 


and 


2.  No  surface  cooling.  Heat  transfer  coefficient,  h = 0. 
For  this  case,  the  boundary  conditions  are: 

f (0,  t)  = £ (a,  t)  = 0 
T (x,  0)  = T 


(9) 


(10) 


which,  with  the  aid  of  Eqs.  (4)  and  (5)  yield: 


oo  -8A2t  cq  e_Wa  2 

T = T + T C e n cos(A  x)  + 0 r'X 

o n n 

n=l 


26t 


K 


f ~ f)sinh(ya) 

SL  SL  Z j 

au 


2 coshpix) 


(11) 


where 

II 

C 

<< 

niT/a 

(12) 

and 

C = 
n 

N ' 
n 

, Eq  e 
-1  0 

K 

■Pa  a 2 2 

— / cos (A  x)  [—  cosh(yx)  - ~ 

0 n y a 

sinh  (ua)]dx  (13a) 

with: 

N = 

a 

/ 

cos^(X 

x)dx  = a/2 

(13b) 

n 

0 

n 

After  integration,  Eq.  (13a)  yields: 

4eq  e"ya  2 , ,,n+l 
C = — ^ ^ — 5—}  sinh(ya) 

n aK  A V+A  2) 

n n 


(13c) 


Note  that  the  boundary  conditions  underlying  Eq.  (11)  as  expressed  by 
Eqs.  (9)  and  (10)  assume  that  the  material  does  not  re-emit  the  absorbed 
radiation.  For  this  reason,  Eq.  (11)  must  be  considered  valid  only  for 
plates  with  low  initial  temperatures  and  over  a time  period  sufficiently 
short  such  that  the  plate  remains  cool  enough  so  that  the  re-emitted 
radiation  to  a first  approximation  can  be  neglected. 

From  Eq.  (11)  it  may  be  noted  that  after  an  initial  transient  the 


5 


temperature  in  the  plate  at  any  point  rises  uniformly  with  time  as 
expressed  by  the  last  term,  with  a temperature  profile  in  the  plate 
given  by  the  third  and  fourth  terms. 


C.  Thermal  Stresses. 

The  general  expression  for  the  thermal  stresses  in  a plate  with  a 
temperature  distribution  which  varies  through-the-thickness  only  and 
which  is  symmetric  about  x = 0,  is^: 

a = [-T  + I Tdx]  (14) 

y,z  ( 1-v)  2a  ' 

where  a is  the  coefficient  of  thermal  expansion,  E is  Young's  modulus, 
and  v is  Poisson's  ratio. 

Substitution  of  Eq.  (7)  or  Eq.  (11)  into  Eq.  (14)  yields  expres- 
sions for  the  thermal  stresses  as  follows: 

1.  Infinite  heat  transfer  coefficient  (h=°°) , or  constant  surface 
temperature. 


4aE(ya)eq  e ya 


. n 


y.z 


00  ~8Xnt.  . (-1)  cos  (A  x) 

o , „ , r (1-e  ) r 1 

cosh(ya)  ~ 5 — t — j T — 


2 2 2 

n=0  (y  +A  ) a A 

n n 


— } (15) 


(l-v)Ka 

Since  ceramic  materials  generally  are  much  weaker  in  tension  than 
in  compression,  the  values  of  the  tensile  stresses  are  of  primary  inter- 
est. Examination  of  Eq.  (15)  indicates  that  the  tensile  stresses  have 
their  maximum  values  at  x = -a,  a and  t = 00  and  can  be  written: 


0 

y.z 


(a,®) 


4aE(ya)eq  e 
o 

(l-v)Ka 


-V  a 


2 


cosh(ya) 


y {(uW)(aA^r1 
n n 

n=0 


(16) 


Eq.  (16)  shows  that  the  magnitude  of  the  thermal  stresses  is  a non- 
linear function  of  the  absorption  coefficient  y. 

Two  limiting  cases  can  be  discussed: 


a.  ya 


For  this  case: 


o (x,°°)  = 0 
y,z 


Eq.  17  implies  that  no  thermal  stresses  are  generated  for  h = °°, 


b.  ya  -*■  0 

For  this  case  Eq.  (16)  can  be  simplified  to: 


0 (a,00)  = 0 (-a,00)  = 

y,Z  y,Z 


64aEy (l-ya)eqQa 
tt4(1-v)K 


using: 


lim  ^ {(y2+X2)(aX2)}_1  = 16a3/TT4  (19) 

ya-K)  n=0  n n 

Eq.  (18)  indicates  that  for  lightly  absorbing  plates  with  ya-K) 
and  h=°°,  the  thermal  stresses  are  a parabolic  function  (with  downward 
curvature)  of  the  absorption  coefficient. 

For  purposes  of  the  selection  of  partially  absorbing  materials  with 
optimum  resistance  to  thermal  fracture  when  subjected  to  radiation 
heating,  Eq.  (18)  for  ya  <<  1 can  be  rewritten  in  terms  of  the  maximum 
radiant  heat  flux,  9inax»  to  which  the  plate  can  be  subjected,  so  that 
the  tensile  stresses  do  not  exceed  the  tensile  strength,  St»  by: 


St(l-v)KTT 

64aEyea2 


(ya  « 1) 


If  the  radiant  heat  flux,  q , is  the  result  of  black-body  radiation 

max 

4 

from  a temperature  source,  T , with  q = pT  , (where  p = Stefan- 

max  max  max 

Boltzmann  constant),  Eq.  (20)  can  be  expressed  also  in  terms  of  the 
maximum  radiation  temperature  to  which  the  plate  can  be  subjected  by: 

T = (— LT)1M  )1/4  (21) 

max  64ap2  “EUe 


Lt. 


2.  No  surface  cooling  (h  = 0) 

Substitution  of  Eq.  (11)  into  Eq.  (14)  yields: 


O = - — l C i 
y,z  1— v , n 
n=l 


eq  aEe  ~ 2 

e nt  cos(Anx)  + Kd'-'v)"  t(f 2 " 27)  sinh  (ya) 

ay 


2cosh(yx) . 

y 


(22) 


The  tensile  thermal  stresses  are  a maximum  at  x = 0 and  t = 00  given 


by: 


„ “ha 

aEeq  e „ _ 

ay,z  = TiVk  {y  + (f  " ^2>slnh^a>} 


(23) 


Two  limiting  cases  of  Eq.  23  can  be  discussed: 

a.  ya  -*■  00 . (24) 

Eq.  (24)  corresponds  to  the  condition  that  all  the  incident  radiant 
heat  is  absorbed  in  the  immediate  surface  layer. 

With  ya  = Eq.  (23)  becomes: 


aEeqQa 


(25) 


y,z  6(l-v)K 

As  expected,  the  thermal  stresses  for  ya  = 00  are  independent  of  the 
absorption  coefficient  and  a function  of  the  emissivity  only, 
b.  ya  « 1. 

no 

For  this  case  Eq.  (23)  with  the  aid  of  the  expansion  e - 1 + ya 
2 2 

+ ^ + ....  can  be  written: 

3 4 

7 aEeq  y (l-ya)a 

°y,z  = 180(1-V)K  (26) 

which  shows  that  for  lightly  absorbing  (uncooled)  plates  with  ya  « 1, 
the  thermal  stresses  are  proportional  to  the  third  power  of  the  absorp- 
tion coefficient. 

Eqs.  (25)  and  (26)  can  be  written  in  terms  of  the  maximum  radiant 


4-  — 
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heat  flux  to  which  the  plate  can  be  submitted,  so  that  the  maximum 
tensile  thermal  stress  does  not  exceed  the  tensile  strength.  In  this 
manner,  for  ya  = 00 , Eq.  (25)  becomes: 


6St(l-v)K 


max  aEea 

which  for  the  maximum  temperature  of  black-body  radiation  becomes: 
6S. (l-v)K 


(27) 


T = { — } 

max  aEeap 


1/4 


(28) 


Similarly  for  ya  « 1,  Eq.  (26)  becomes: 


180St(l-v)K 
?Tnax  = 7aEey3a4 


or  in  terms  of  T , Eq.  (29)  becomes: 
max 


(29) 


T = {• 
max 


180St(l-v)K 

“ TT 

7aEepy  a 


-} 


1/4 


(30) 


III.  NUMERICAL  RESULTS,  DISCUSSION  AND  CONCLUSIONS 

Figure  1 shows  the  maximum  value  of  the  non-dimensional  tensile 
thermal  stress  as  a function  of  the  product  ya  for  the  two  values  of 
heat  transfer  coefficient  h = 0,  °°,  considered  in  the  analysis.  For 
both  cases  (h  = 0,  00 ) the  thermal  stresses  are  identically  equal  to 
zero,  for  ya  = 0.  For  h = 0,  the  value  of  maximum  tensile  thermal 
stress  increases  with  ya,  with  a maximum  value  of  0.167,  as  ya  °°.  In 
contrast,  for  h = 00 , the  values  of  maximum  non-dimensional  stress  go 
through  a maximum  at  ya  *=  1.344  at  a value  of  non-dimensional  stress 
equal  to  0.279  followed  by  a decrease  to  zero  as  ya  -*■  00 . The  latter 
effect  arises  because  as  ya  -*■  00  the  absorbed  radiant  energy  is  trans- 
formed into  heat  in  the  immediate  surface  regions  of  the  plate  to  be 
removed  immediately  by  the  surrounding  medium.  For  this  reason  as 
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ya  -*■  00  no  temperature  increase  in  the  plate  will  occur  and  the  plate 


will  remain  stress-free. 

It  is  critical  to  note  that  the  present  analysis  was  based  on  the 
assumption  that  the  temperature  of  the  plate  when  the  stresses  have 
reached  a value  near  their  final  value  (Fig.  1)  is  still  sufficiently 
low  that  the  re-emitted  radiation  to  a first  approximation  can  be  neg- 
lected. The  validity  of  this  assumption  was  established  by  calculating 
the  transient  temperature  and  thermal  stresses  as  a function  of  time  for 
a range  of  values  of  pa.  For  h = 00 , the  values  for  the  temperatures  and 
stresses  as  a function  of  time  are  shown  in  Figs.  2 and  3,  and  for  h = 0 

in  Figs.  4 and  5,  respectively.  For  h = 0,  the  temperatures  are  unbounded 

8 

as  t + ®.  For  this  reason,  as  done  in  a previous  study  , the  tempera- 
tures in  the  plate  were  calculated  for  the  time  at  which  the  stresses 

have  reached  95%  of  their  maximum  value  corresponding  to  q for  a 

max 

given  numerical  example.  For  values  of  h = 00  and  Via  = 1.344  and  h = 0 
and  ya  “ 2.5,  these  temperatures  together  with  the  physical  property 
values  chosen  are  listed  in  Table  I.  These  property  values  correspond 
approximately  to  those  of  a high-density  polycrystalline  alumina.  The 
numerical  results  indicate  clearly  that  the  re-emitted  radiation  is  a 
very  small  fraction  of  the  incident  radiation.  Similar  quantitative 
results  can  be  obtained  for  other  values  of  ya. 

For  h = °°,  the  calculated  value  of  the  temperature  in  the  plate  at 

2 

3t/a  ■ 1 also  is  close  to  the  final  steady-state  temperature  attained 
2 

as  3t/a  -*  °°,  so  that,  at  least  for  the  numerical  example  given  in  Table 
I,  the  radiative  heat  losses  need  not  be  considered  at  all.  For  h = 0, 
however,  after  the  initial  transient  the  temperatures  will  rise  linearly 
in  time,  so  that  after  a sufficient  time  period  has  elapsed,  radiation 
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heat  losses,  thermo-visco-elastic  stress  relaxation  or  melting  phenomena 
need  be  taken  into  account.  However,  as  indicated  by  the  numerical 
example  given  in  Table  I,  if  thermal  stress  fracture  occurs  it  will  take 
place  during  the  initial  transient  before  the  plate  is  hot  enough  to  re- 
emit the  absorbed  heat  flux  and  at  temperatures  at  which  the  plate  is 
still  brittle. 

In  general,  however,  it  should  be  noted  that  the  validity  of  neg- 
lecting the  re-emitted  radiation  should  be  examined  for  each  individual 
case  under  consideration.  For  instance,  for  ultra-strong  structures  or 
those  with  low  values  of  coefficient  of  thermal  expansion,  the  heat 
fluxes  required  for  failure  may  be  so  high  with  corresponding  high  plate 
temperatures  that  the  re-emitted  radiation  may  be  an  appreciable  frac- 
tion of  the  incident  heat  flux.  If  so,  the  expressions  for  the  stresses 
represent  overestimates,  such  that  the  equations  for  q or  T are 
conservative.  In  principle,  solutions  for  the  temperatures  and  stresses 
for  high  levels  of  re-emitted  radiation  can  be  obtained  by  numerical 
methods.  Such  numerical  calculations  also  may  include  the  effect  of 
temperature  on  the  pertinent  material  properties,  multiple  reflections, 
the  spectral  dependence  of  the  emissivity  and  absorption  coefficient,  as 

well  as  non-linear  effects  such  as  thermo-visco-elastic  stress  relaxa- 

12  13 

tion.  Such  numerical  evaluations  as  judged  by  other  studies  ’ are 
more  complex  than  the  present  analysis.  The  advantage  of  the  present 
analysis  in  spite  of  its  simplifying  assumptions  is  that  the  stresses 
are  expressed  analytically  in  terms  of  the  appropriate  variables,  which 
permits  a direct  assessment  of  the  relative  role  of  each  material  property 
in  establishing  the  magnitude  of  thermal  stress. 

In  this  respect,  it  is  of  interest  to  note  the  role  of  absorption 
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in  determining  the  magnitude  of  thermal  stress.  For  ya  « 1 and  h = 00 
the  thermal  stresses  are  linearly  proportional  to  the  absorption  coef- 
ficient. In  contrast,  for  pa  <<  1 and  h = 0,  the  thermal  stresses  are 
proportional  to  the  cube  of  the  absorption  coefficient.  Coupled  with 
the  role  of  the  absorption  in  governing  the  magnitude  of  thermal  stress 
is  also  an  unusual  size  effect.  For  h = 00  and  ya  « 1,  the  thermal 
stresses  are  proportional  to  the  square  of  the  plate  thickness.  For  h = 
0 and  ya  <<  1,  the  maximum  thermal  stresses  are  proportional  to  the 

fourth  power  of  the  plate  thickness.  These  size  effects  are  in  contrast 

3 

to  those  for  convective  heat  transfer  , for  which  at  low  Biot  number 
(ah/K  < 1)  the  thermal  stresses  are  directly  proportional  to  the  size 
and  independent  of  size  for  high  Biot  number  (ah/K  -*■  °°) . 

The  expressions  developed  for  the  maximum  permissible  incident  heat 
flux  and  black-body  radiation  permit  the  formulation  of  additional 
"thermal  stress  resistance"  parameters,  which  can  be  used  for  the  selec- 
tion of  the  material  with  the  optimum  combination  of  properties  for 
maximum  thermal  stress  resistance.  These  parameters  can  be  defined  as 
follows: 

For  h = 00  and  ya  « 1,  Eqs.  (20)  and  (21)  yield: 


S (1-V)K  S (1-V)K  1/4 

and  {"  oEe'y"  } 


Similarly,  for  h = 0 and  ya  « 1,  Eqs.  (29)  and  (30)  yield: 

S (1-V)K  S (l-v)K  1/4 

{— —)  and  {-£ r} 

aEeyJ  aEeyJ 

Furthermore,  for  h = 0,  ya  -*■  00 , Eq.  (27)  yields: 

S „(1-V)K 


■ — - jL  ..q 


The  parameter  of  Eq.  (33)  raised  to  the  1/4  power  was  presented  in 

g 

an  earlier  study  . 

With  these  five  additional  parameters,  the  number  of  non-redundant 
thermal  stress  resistance  parameters  available  for  material  selection 
now  has  grown  to  a total  near  thirty.  Also,  with  the  introduction  of 
the  absorption  coefficient  the  number  of  material  properties  presently 
known  to  affect  the  thermal  stress  resistance  of  brittle  materials  has 
reached  a total  of  eighteen. 

From  the  point  of  view  of  material  selection  in  order  to  minimize 
the  incidence  of  thermal  stress  fracture,  it  should  be  noted  that  for  h 
= 0,  a partially  transparent  ceramic  is  expected  to  have  greater  thermal 
stress  resistance  than  an  opaque  ceramic.  On  the  other  hand,  for  h = 00 , 
low  thermal  stresses  can  be  achieved  for  either  highly  transparent 
ceramics  (pa  0)  or  for  opaque  materials  with  ya  -*  °°. 

In  summary,  an  analysis  has  been  presented  for  the  thermal  stress 
resistance  of  partially  absorbing  brittle  ceramics  in  the  form  of  flat 
plates  subjected  to  symmetric  normally  incident  radiation  heating  for 
two  specific  convective  cooling  conditions.  Appropriate  thermal  stress 
resistance  parameters  were  derived. 
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TABLE  I.  Calculated  Maximum  Heat  Flux  and  Temperature 


(ABSORBTION  COEFFICIENT) (PLATE  HALF  THICKNESS) 
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Fig.  4.  Transient  temperatures  at  surface  of  partially  absorbing  flat  plate  symmetrically  heated 
by  normally  incident  radiation  and  cooled  by  convection  with  heat  transfer  coefficient, 


0.20 


Fig.  5.  Transient  tensile  thermal  stresses  at  center  of  partially 
absorbing  flat  plate  symmetrically  heated  by  normally  in- 
cident radiation  and  cooled  by  convection  with  heat  trans- 
fer coefficient,  h = o. 


CHAPTER  II. 

ANALYSIS  OF  THERMAL  STRESS  RESISTANCE 
OF  MICRO-CRACKED  BRITTLE  CERAMICS. 


— - - k 


ANALYSIS  OF  THERMAL  STRESS 


RESISTANCE  OF  MICRO-CRACKED  BRITTLE  CERAMICS 


D.  P.  H.  Hasselman  and  J.  P.  Singh 
Department  of  Materials  Engineering 
Virginia  Polytechnic  Institute  and  State  University 
Blacksburg,  Virginia  24061 


' 


f 


ABSTRACT 

An  analysis  is  presented  for  the  effect  of  micro-cracking  on  thermal 
stress  resistance  of  brittle  ceramics.  Expressions  are  obtained  or  de- 
rived for  thermal  conductivity,  elastic  properties,  coefficient  of  thermal 
expansion  and  strength  in  terms  of  crack-size,  crack  density  and  crack 
interaction  for  a micro-cracked  material.  It  is  shown  that  the  values  of 
Young  modulus,  Poisson's  ratio,  coefficient  of  thermal  expansion  and 
thermal  conductivity  decrease  with  increasing  degree  of  micro-cracking. 

Strength,  however,  may  increase  or  decrease  depending  upon  the  nature  of 

: 

crack  interaction.  The  combined  effect  of  micro-cracking  on  the  above 
mentioned  properties  results  in  a significant  increase  in  thermal  stress 
resistance.  In  particular,  it  makes  a micro-cracked  material  very  attractive 
for  the  application  where  high  thermal  stress  resistance  with  good  thermal 

» 

insulating  property  is  required. 


I. 


Introduction 


Brittle  ceramic  materials  for  structural  applications  at  high  tem- 
perature generally  are  highly  susceptible  to  catastrophic  failure  under 

conditions  of  thermal  stress.  As  demonstrated  experimentally^  and  analy- 
2 3 

tically  ’ , materials  with  high  densities  of  micro-cracks  may  offer  a 
solution  to  the  low  thermal  stress  resistance  of  brittle  materials.  Micro- 

cracked  materials  exhibit  high  strains-at-f racture  high  fracture 

6 2 
toughness  and  on  thermal  failure  should  undergo  stable  crack  propagation  ' 

rather  than  unstable  (catastrophic)  crack  propagation  encountered  for  many 

homogeneous  high-strength  ceramics.  The  high  thermal  stress  resistance 

of  micro-cracked  materials  primarily  can  be  attribured  to  the  effect  of 

the  micro-cracks  on  elastic  behavior.  Heavily  micro-cracked  materials, 

in  fact,  can  exhibit  elastic  moduli  substantially  lower  than  the  non- 

4 5 

micro-cracked  materials  ’ . Unfortunately  micro-cracking  can  cause  a 

significant  decrease  in  thermal  conductivity,  with  a corresponding  de- 

7,8,9 

crease  in  thermal  stress  resistance 

A quantitative  assessment  of  the  effect  of  micro-cracking  on  thermal 
stress  resistance  requires  an  analysis  of  the  combined  effects  of  micro- 
cracking on  elastic  behavior,  thermal  conductivity,  coefficient  of 
thermal  expansion  and  the  fracture  stress.  This  latter  property,  in 
particular,  may  be  strongly  affected  by  crack  interaction. 

It  is  the  purpose  of  this  study  to  present  literature  information 
on  the  effect  of  micro-cracking  on  the  principal  properties  which  af- 
fect thermal  stress  resistance  or  to  derive  such  information,  if  required. 
From  the  combined  results,  an  assessment  then  is  made  of  the  effect  of 
micro-cracking  on  thermal  stress  resistance. 


I 

i 
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II.  Analysis 

A.  General  Approach  and  Boundary  Conditions. 

For  the  information  on  the  effect  of  micro-cracks  on  material 
properties,  only  those  expressions  will  be  presented  or  derived  which 
are  sufficient  for  an  analysis  of  the  effect  of  micro-cracking  on  ther- 
mal stress  resistance.  For  a review  and  critical  analysis  of  the 
mathematical  approaches  to  and  the  assumptions  underlying  the  equations 
given,  the  reader  is  referred  to  the  cited  literature. 

As  will  be  shown,  the  relative  effect  of  micro-cracking  on  a given 
material  property,  q is  derived  in  the  literature  or  in  the  present 
paper  in  the  general  form:  q = qQ(l-£)  or  q = q^Cl+e)  \ where  qQ  is 
the  value  of  the  material  property  of  the  crack-free  material.  The 
quantity  e,  which  includes  a constant  is  a function  of  the  crack  den- 
sity and  crack  size.  Since  for  e <<  1,  (1-e)  = (1+e)  \ for  simplicity 
of  comparison,  all  equations,  whenever  appropriate  will  be  given  in  the 
form  of  q = qo(l+e)  \ This  also  avoids  the  problem  of  obtaining  zero 
or  negative  property  values  for  e - 1. 

In  order  to  reduce  the  total  number  of  solutions  which  could  be 
presented  to  a manageable  total,  for  so  far  as  possible  the  analysis  will 
be  limited  to  a crack  geometry  in  the  form  of  ellipsoids  of  revolution 
of  very  low  aspect  ratio  (i.e.,  cracks  which  in  the  limit  approach  zero 
thickness).  Also,  all  cracks  will  be  assumed  to  be  of  uniform  size  and 
geometry. 

Finally,  also  it  will  be  assumed  that  the  crack  sizes  and  tempera- 
ture gradients  involved  during  the  thermal  stress  are  sufficiently  small, 
so  that  crack  instability  due  to  the  temperature  non-uniformity  in  the 
vicinity  of  the  cracks  does  not  occur.  As  is  easily  ascertained  from  the 


theory  of  Goodier  and  Florence  , this  latter  assumption  is  easily 
satisfied  for  cracks  of  the  order  of  Griffith  flaws  and  temperature 
gradients  of  the  order  of  those  usually  found  in  ceramic  technology. 


B.  Effect  of  Micro-cracks  on  Thermal  Conductivity 

Solutions  for  the  effect  of  micro-cracks  in  the  shape  of  ellipsoids 

11,12 

of  revolution  of  zero  thickness  on  thermal  conductivity  are  available 
as  follows: 

1.  Micro-cracks  with  statistical  random  orientation. 

The  effect  of  micro-cracks  with  statistical  random  orientation  on 

12 

thermal  conductivity  can  be  expressed  : 


K = K (1  + 8Nb3/9)_1 
o 


where  K,  K are  the  thermal  conductivities  of  the  micro-cracked  and  crack- 
o 

free  material,  resp.,  N is  the  number  of  cracks  per  unit  volume  and  b 

is  the  radius  of  revolution  of  the  cracks. 

2.  Micro-cracks  oriented  with  plane  of  crack  propagation 
perpendicular  to  the  direction  of  heat  flow. 

The  effect  of  micro-cracks  with  this  preferred  orientation  can  be 


......  11,12 

derived  to  be  : 


K = K (1  + 8Nb3/3)  1 
o 


3.  Micro-cracks  oriented  with  plane  of  cracks  parallel  to 
heat  flow. 

12 

For  this  crack  orientation  : 

K - K 


The  result  of  eq.  3 is  reasonable  since  the  cracks  are  not  expected 
to  impede  the  heat  flow.  Figure  1 compares  the  relative  thermal  conductivity 


for  the  three  orientations  as  a function  of  the  quantity  Nb  . Cracks 


L 


•» Lx 


with  preferred  orientation  perpendicular  to  the  direction  of  heat  flow 
have  the  highest  relative  effect  on  thermal  conductivity,  followed  by 
randomly  oriented  cracks  and  cracks  parallel  to  the  direction  of  heat 

flow.  Experimental  data  for  the  effect  of  micro-cracking  on  thermal 

7 8 3 3 

conductivity  ’ suggest  that  Nb  = 2,  so  that  the  range  of  0 < Nb  <5 

chosen  for  figure  1 must  be  considered  realistic. 
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C.  Effect  of  Cracks  on  Elastic  Properties 

This  subject  has  received  considerable  theoretical  attention-1 
including  the  effect  of  the  absence  and  the  presence  of  crack  interaction. 
As  for  thermal  conductivity,  crack  orientation  plays  a significant  role 
in  the  effect  of  micro-cracks  on  elastic  behavior,  as  presented  or  derived 
as  follows: 

1.  Random  orientation 

For  the  purpose  of  the  present  study,  convenient  expressions  for  the 
effect  of  randomly  oriented  micro-cracks  on  Young's  modulus  and  Poisson's 
ratio  were  obtained  by  Salganik-^  as  follows: 


16(10-3v  )(l-v  ) ~ . 

E = E [1  + , ° Nb  J 

o 45(2-v  ) 


(4) 


16(3-v  ) (l-v  ) 3 

» ' “o'1  + - "IS (2-0  — " I 


(5) 


where  E and  v are  Young's  modulus  and  Poisson's  ratio  of  the  micro-cracked 

material,  and  E and  v are  Young's  modulus  and  Poisson's  ratio  of  the 
o o 

crack-free  material. 

3 

Figures  2 and  3 show  the  dependence  of  E/Eq  and  v on  the  quantity  Nb  . 
Note  that  to  a first  approximation,  E/Eq  is  independent  of  the  value  of 
vq  and  can  be  expressed: 


4 


(6) 


E = E [1  + 16Nb3/9]~1 
o 


For  purpose  of  the  discussion  of  the  effect  of  micro-cracking  on 


the  coefficient  of  thermal  expansion,  an  expression  for  the  effect  of 

14 

micro-cracking  on  bulk  modulus  is  required.  As  derived  by  Walsh 
(also  directly  obtainable  from  eqs.  4 and  5)  this  effect  can  be  express- 


16(l-v  2) Nb2  , 

B - V1  * 9(1-2^) r <7) 

2.  Micro-cracks  with  crack  plane  oriented  perpendicular  to 
uni-axial  stress  direction. 

For  circular  micro-cracks  oriented  perpendicular  to  the  stress  direc- 
tion, Young's  modulus  in  the  direction  of  stress  can  be  derived  from  the 
20 

solution  of  Sack  which  shows  that  for  a single  crack  the  change  in 
potential  energy  (AW)  is: 

AW  = 8(l-v  2)o2b3/3E  (8) 

o o 

where  E is  Young's  modulus  of  the  crack-free  material.  For  N cracks  per 
o 

unit  volume,  assuming  the  absence  of  crack  interaction,  the  change  in  poten- 
tial energy  is: 

AW.,  = 8(l-v  2)o2Nb3/3E  (9) 

No  o 

Since  the  potential  energy  per  unit  volume  in  the  crack-free  material 

equals  a / 2E  , the  total  potential  energy  in  the  material  with  cracks: 
o 


W = o2/2E  +8(l-v  2)o2Nb3/3E 
too  o 


In  terms  of  the  effective  Young's  modulus  (E)  of  the  material  the  poten- 


tial energy  (Wf)  becomes: 


Wt  * o /2E 


r 
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Equating  eqs.  10  and  11  yields: 

E = E [1+16 (1-v  2) Nb 3/ 3 ] _ 1 (12) 

o o 

3.  Micro-cracks  with  crack  plane  parallel  to  stress  direction. 

For  this  case: 

E = E (13) 

o 

21 

as  can  be  derived  easily  from  the  upper  bound  for  the  Young's  modulus 
of  a composite  with  a dispersed  phase  with  Young's  modulus  and  volume 
fraction  which  approach  zero. 

The  above  expressions  for  the  effect  of  micro-cracking  on  elastic 
behavior  strictly  are  valid  for  dilute  concentrations  of  cracks  for  which 
crack  interaction  effects  can  be  neglected.  At  higher  concentrations 
at  which  the  distance  between  cracks  is  on  the  order  of  or  less  than  the 
crack  size,  crack  interaction  effects  undoubtedly  can  play  a major  part. 
Unfortunately,  as  far  as  the  present  writers  are  aware,  no  solutions 
for  interaction  effects  between  cracks  in  the  form  of  ellipsoids  of  re- 
volution on  elastic  behavior  have  appeared  in  the  literature.  Nevertheless, 

3 15  18 

such  solutions  are  available  for  arrays  of  cracks  in  plates  * ’ , which 

can  give  a qualitative  or  semi-quantitative  indication  of  the  interaction 
effects  between  ellipsoidal  cracks. 

22 

From  the  solution  of  Yokobori  and  Ichikawa  for  co-planar  rows  of 
cracks,  Hasselman*5  derived  the  elastic  behavior  of  a rectangular  array 
of  cracks  depicted  in  figure  4,  in  which  interaction  effect  between  the 
rows  was  assumed  absent.  The  effective  Young's  modulus  (E)  of  the  plate 
perpendicular  to  the  plane  of  cracks  (for  plane  strain)  is: 

16Nd2(l-v  2) 

E = E [1 — Hn  cos  rrr  ] (14) 

o it  2d 
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where  a is  the  half-length  of  the  cracks,  N is  the  density  of  the  cracks 
per  unit  area  and  d is  defined  in  figure  4. 

For  very  dilute  concentrations  of  cracks  (N  o,  d ->•  00 ) eq.  14 
becomes: 

E = E [1  + 2ttN(1-v  2)a2]_1  (15) 

o o 

3 

which  corresponds  to  the  solution  obtained  by  Hasselman  for  dilute  con- 
centration (absence  of  crack  interactions)  adjusted  for  the  plane  stress 
condition. 

Comparison  between  eqs.  14  and  15  indicates  that  interaction  between 
the  co-planar  cracks  causes  Young's  modulus  to  drop  to  a value  below 
the  corresponding  value  for  lack  of  crack  interaction. 

The  effect  of  interaction  between  the  parallel  rows  of  the  co-planar 

cracks  on  Young's  modulus  can  be  ascertained  from  the  results  of  Delameter, 
18 

Herrman  and  Barnett  , shown  in  figure  5,  which  combine  the  effects  of 
interaction  between  the  rows  of  co-planar  cracks  and  between  the  co-planar 
cracks.  As  expected, for  a given  value  of  2a/c,  Young's  modulus  decreases 
with  increasing  2a/d.  However,  for  a given  value  of  2a/d,  the  uoward 
curvature  of  E/EQ  with  increasing  2a/c,  shows  that  the  nature  of  interaction 
between  the  rows  of  co-planar  cracks  is  such  that  Young's  modulus  is 
increased;  in  other  words,  the  interaction  effect  between  neighboring 
rows  of  cracks  is  to  reduce  the  effect  on  Young's  modulus  due  to  an  in- 
dividual non-interacting  crack.  This  observation  will  be  shown  to  be 
critical  to  the  effect  of  the  micro-cracks  on  strength. 


D.  Strength 

The  introduction  of  micro-cracking  in  a material  is  expected  to  have 
an  effect  on  strength.  The  micro-cracks  themselves  can  constitute  the 
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failure-initiating  flaws  or  may  have  ar.  effect  on  pre-existing  larger 
flaws. 


I 


I 

I 

1 

X 

X 

l 

I 

I 

I 


For  dilute  concentrations  of  micro-cracks  for  which  crack  inter- 
actions effects  are  absent,  the  strength  will  be  independent  of  the 

20 

crack  density.  As  derived  by  Sack  for  a circular  crack  the  critical 
fracture  stress  (St>  is: 

St  = [TrEYf/2a(l-vo2)]1/2  (16) 

in  which  y^  is  the  surface  fracture  energy  per  unit  area  of  the  newly 
created  fracture  surface. 

At  high  crack  densities,  however,  crack  interaction  effects  on 

strength  cannot  be  ignored.  To  examine  this  effect,  reliance  must  be 

placed  on  solutions  for  cracks  in  plates.  For  a single  co-planar  row 
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of  cracks,  the  fracture  stress  as  derived  by  Yokobori  and  Ichikawa  can 
be  expressed  (for  conditions  of  plane  strain)  as: 


ey. 


St  = [ 


2(1-vq  )d 


„ ,ira.1l/2 
c°t  (^)] 


(17) 


Eq.  17,  as  d + «,  approaches  the  Griffith  solution  for  a single 


crack: 


Sc  = [EYf/*a(l-vo2)  ]1//2 


(18) 


Comparison  of  eq.  17  and  18  shows  that  the  interaction  between  co- 

planar  cracks  is  to  cause  a decrease  in  the  value  of  strength  from  the 

value  corresponding  to  a single  crack. 

The  effect  of  crack  interactions  between  rows  of  co-planar  cracks  (as 

18 

shown  in  figure  4)  was  calculated  by  Delameter,  Herrman  and  Barnett 
Their  results  in  terms  of  relative  stress  intensity  factor  are  shown  in 
figure  6.  These  results  show  that  for  a given  value  of  the  distance  (d) 


i 
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between  co-planar  cracks,  the  stress  intensity  factor  decreases  with  de- 
creasing distance  between  rows  of  co-planar  cracks,  which  should  lead  to 
an  increase  in  strength.  This  leads  to  the  surprising  result  that  a 
material  with  many  cracks  of  the  proper  relative  orientation  and  inter- 
crack spacing  can  be  stronger  than  the  non-microcracked  material. 

On  the  other  hand,  as  indicated  by  eq . 17  and  the  results  shown  in 
figure  6 for  small  values  of  d/2a,  the  interaction  between  cracks  can  be 
such  that  strength  is  decreased.  In  part,  this  explanation  serves  to 
explain  the  generally  lower  strength  of  micro-cracked  materials  in  com- 
parison to  their  non-microcracked  counterparts. 

Another  explanation  for  a decrease  in  strength  can  be  based  on  the 
fact  that  failure  of  brittle  materials  can  be  initiated  from  pre-existing 
flaws  of  sizes  much  larger  than  the  micro-cracks.  In  this  case,  some 
micro-cracks  which  form  may  lie  co-planar  with  the  pre-existing  flaw, 
which  in  effect  experiences  an  increase  in  its  effective  length  and 
causes  a corresponding  decrease  in  strength.  In  fact,  it  is  thought  that 
this  latter  mechanism  is  primarily  responsible  for  the  decrease  in  strength 
of  brittle  materials  which  exhibit  extensive  micro-cracking. 

E.  Effect  of  Cracks  on  Coefficient  of  Thermal  Expansion 

The  coefficient  of  thermal  expansion  of  a homogeneous  material  does 
not  depend  on  the  geometry  of  the  material.  For  this  reason,  it  is  not 
expected  that  a brittle  material  with  cracks  of  zero  width  will  have  a 
coefficient  of  thermal  expansion  different  from  the  non-cracked  material. 
For  materials  which  are  heterogeneous  in  thermal  expansion  behavior,  the 
above  conclusion  is  not  expected  to  be  valid.  In  such  materials,  such  as 
composites  with  a mismatch  between  the  coefficients  of  thermal  expansion 
of  the  individual  components  or  in  polycrystalline  materials  which  exhibit 
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thermal  expansion  anisotropy,  internal  stresses  can  arise  which  can 
control  the  magnitude  of  the  overall  coefficient  of  thermal  expansion 
of  the  composite  or  polycrystalline  aggregate. 

In  general,  the  coefficient  of  thermal  expansion  of  such  hetero- 
geneous materials  is  a function  of  the  coefficients  of  thermal  expansion, 
the  volume  fraction  and  the  elastic  moduli  of  each  component.  Micro- 
cracking as  the  result  of  the  internal  stresses  will  change  the  elastic 
moduli  with  a resulting  change  in  the  overall  coefficient  of  thermal 
expansion.  Quantitatively,  this  effect  will  depend  on  the  distribution 
of  each  component,  the  internal  stress  with  the  heterogeneous  system 
and  other  microstructural  variables.  Some  general  features,  however, 
of  the  effect  of  micro-cracking  on  the  coefficient  of  thermal  expansion 
of  heterogeneous  materials  can  be  examined  on  the  basis  of  an  expression 

for  the  coefficient  of  thermal  expansion  of  a composite  (ac)  derived  by 
23 

Turner  , which  for  a two-component  system  can  be  written: 


a 

c 


alVlBl  + a2V2B2 
V1B1  + V2B2 


(19) 


where  a is  the  coefficient  of  thermal  expansion,  B is  the  bulk  modulus 
and  V is  the  volume  fraction  of  each  component  1 and  2.  Generally,  for 
materials  manufactured  at  high  temperature,  on  cooling  the  component  with 
the  higher  coefficient  of  thermal  expansion  will  be  subjected  to  a high 
level  of  residual  tensile  stress.  These  stresses  when  of  sufficient 
magnitude  can  cause  extensive  micro-cracking  and  a corresponding  decrease 
in  bulk  modulus.  If  so,  examination  of  eq.  19  reveals  that  in  a composite 
in  which  one  of  the  components  develops  micro-cracks,  the  overall  coef- 
ficient of  thermal  expansion  (ac)  will  approach  the  value  for  the  component 


10 


with  the  lower  coefficient  of  thermal  expansion.  Similar  considerations 
apply  to  the  coefficient  of  thermal  expansion  of  a polycrystalline  ag- 
gregate of  a material  which  exhibits  anisotropy  in  thermal  expansion  be- 
havior. This  explains  the  very  low  or  even  negative  values  for  the 
coefficient  of  thermal  expansion  of  micro-cracked  polycrystalline  materials, 
in  spite  of  the  fact  that  the  value  of  the  coefficient  of  thermal  expansion 
averaged  over  crystallographic  orientation  still  can  be  quite  appreciable. 
Figure  7 shows  the  effect  of  micro-cracking  on  the  relative  coefficient 
of  thermal  expansion  of  a hypothetical  composite  with  B^  = B^,  = 

0.5  and  <*2  = 4(1^.  It  can  be  noted  that  the  micro-cracking  for  reasonable 
3 

values  for  Nb  can  reduce  the  coefficient  of  thermal  expansion  to  half  the 

value  of  the  non-microcracked  material.  The  above  conclusions  are  in 
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general  agreement  with  the  observations  of  Rossi 

F.  Effect  of  Micro-Cracking  on  Thermal  Stress  Resistance 
The  relative  effect  of  micro-cracking  on  thermal  stress  resistance 
by  the  simultaneous  effects  of  the  micro-cracking  on  the  thermal  conducti- 
vity, Young's  modulus,  Poisson's  ratio,  strength  and  the  coefficient  of 

thermal  expansion,  can  be  judged  most  conveniently  on  the  basis  of  the 
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well-known  thermal  stress  resistance  parameters  : 

R = St(l-v)/aE,  R'  = St(l-v)K/aE  (20) 

For  heterogeneous  materials  such  as  composites  or  micro-cracked 
materials  with  a micro-structural  scale  much  smaller  than  the  body  size, 
the  individual  properties  involved  in  R and  R'  correspond  to  those  of  the 
composite  or  micro-cracked  material  as  a whole. 

Materials  with  high  values  of  the  parameters  R and  R'  should  exhibit 
good  thermal  stress  resistance.  Accordingly,  for  the  present  case,  if 


micro-cracking  increases  thermal  stress  resistance,  the  combined  effect 
of  the  micro-cracking  on  the  individual  properties  should  be  such  that 
an  increase  in  the  parameters  R and  R'  is  achieved. 

3 

Figure  8 shows  the  variation  of  R and  R*  with  Nb  for  randomly 
oriented  cracks  calculated  from  the  combined  effects  of  micro-cracking 
on  thermal  conductivity.  Young's  modulus,  Poisson's  ratio  and  the  coef- 
ficient of  thermal  expansion  as  shown  in  figures  1,  2,  3 and  7 with  the 
assumption  that  crack-interaction  effects  are  such  that  the  strength 
remains  invariant.  For  crack-interaction  effects  such  that  the  stress 
intensity  factor  decreases  (i.e.,  strength  increases)  as  shown  in 
figure  6,  the  corresponding  values  for  R will  lie  above  the  ones  shown 
in  figure  8.  Even  for  crack  interactions  such  that  strength  is  decreased, 
an  increase  in  R is  expected  since  the  combined  effect  of  micro-cracking 
on  strength  and  Young's  modulus  at  least  for  moderate  values  of  crack 

density  is  such  as  to  increase^  the  ratio  of  S^/E.  This  latter  observa- 

4 5 

tion  is  well  documented  by  experimental  ’ observations.  An  increase 
in  St/E  will  lead  to  an  increase  in  R irrespective  of  additional  changes 
in  Poisson's  ratio  and  the  coefficient  of  thermal  expansion.  Since 
crack  interactions  also  are  expected  to  affect  the  thermal  conductivity, 
no  unambiguous  statement  with  respect  to  the  changes  in  R'  can  be  made, 
until  theoretical  solutions  for  this  effect  become  available.  However, 
since  Young's  modulus  is  affected  more  strongly  by  micro-cracking  than  ther- 
mal conductivity,  it  is  reasonable  to  expect  that  with  crack-interaction 
effects  taken  into  account,  the  value  of  R'  will  not  drop  below  the  value 
for  the  non-cracked  material.  Regardless  of  the  details  of  the  crack 
geometry  and  interactions,  the  present  analytical  results  clearly  sub- 
stantiate the  experimental  observation^  that  randomly  oriented  micro-cracks 


can  lead  to  significant  increases  in  thermal  stress  resistance. 
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Significant  increases  in  thermal  stress  resistance  should  be  ob- 
tainable, at  least  in  principle  for  micro-cracks  with  preferred  orientation. 
Generally,  in  cases  involving  severe  heat  flux,  at  least  near  surfaces 
the  maximum  values  of  thermal  stress  arise  in  a direction  perpendicular 
to  the  direction  of  heat  flow.  For  this  reason,  it  is  desirable  to 
maintain  a value  of  thermal  conductivity  in  the  direction  of  heat  flow 
as  high  as  possible.  At  the  same  time.  Young's  modulus  and  the  coeffi- 
cient of  thermal  expansion  perpendicular  to  the  direction  of  heat  flow 
should  be  kept  as  low  as  possible.  This  condition  can  be  achieved  by  pro- 
moting oriented  micro-cracking  such  that  the  plane  of  the  micro-cracks 
lies  parallel  to  the  direction  of  heat  flow,  with  random  or  preferred 
orientation  perpendicular  to  the  heat  flow.  In  this  manner,  the  effective 
thermal  conductivity  in  the  direction  of  heat  flow  as  indicated  by  eq.  3 
is  unaffected  whereas  Young's  modulus  (eq.  12)  and  the  coefficient  of 
thermal  expansion  in  the  direction  of  maximum  thermal  stress  can  be  re- 
duced. 

A non-uniform  distribution  of  cracks  of  random  or  preferred  orienta- 
tion can  also  be  used  to  advantage  in  improving  thermal  stress  resistance. 
Such  a non-uniform  distribution  of  cracks  results  in  a spatial  variation 
in  thermal  conductivity,  which  under  the  proper  conditions  will  lead  to 

a decrease  in  maximum  value  of  tensile  thermal  stress  and  a corresponding 

26 

increase  in  thermal  stress  resistance 

The  effect  of  micro-cracking  on  thermal  stress  resistance  can  be 
regarded  as  beneficial  from  another  viewpoint  as  well.  High-temperature 
materials  often  are  expected  to  perform  as  thermal  insulators  for  the 
purpose  of  the  containment  of  high  temperatures.  Such  materials  then 
may  have  the  requirement  of  high  thermal  stress  resistance  in  combination 

13 


with  low  thermal  conductivity.  This  latter  requirement  is  incompatible 
with  the  first  which  generally  requires  high  thermal  conductivity.  This 
incompatibility  may  be  satisfied  by  using  heavily  micro-cracked  materials. 

The  requirement  of  high  ratio  of  thermal  stress  resistance  to  thermal 
conductivity  can  be  expressed  in  terms  of  "f igures-of-raerit"  which  are 
derived  by  dividing  the  thermal  stress  resistance  parameters  R and  R' 

(eq.  20)  by  the  thermal  conductivity  which  yields: 

Rr  = St(l-v)/aEK  and  = St(l-v)/aE 

Note  that  the  f igure-of-merit  R^  dimensionally  is  equivalant  to  the  para- 
meter R.  It  is  demonstrated  earlier  that  due  to  their  effect  on  Young's 
modulus,  micro-cracks  can  cause  a major  increase  in  the  parameter  R,  and 
therefore  also  in  the  f igure-of-merit  R^.  This  effect  is  even  more  pro- 
nounced for  the  f igure-of-merit  R , which  as  the  result  of  the  simultaneous 
decrease  in  E and  K,  should  show  a major  increase  as  the  result  of  the 
micro-cracking . 

For  the  above  derivations  and  discussions  the  crack  density  and 
crack  size  were  assumed  to  be  constant.  However,  in  micro-cracked 
materials  the  number  and  the  size  of  the  cracks  is  a function  of  the 
degree  of  cooling.  Furthermore,  due  to  crack-closure  and  crack  healing, 
reheating  of  micro-cracked  materials  can  result  in  the  recovery  of  many 
of  the  physical  properties  close  to  the  non-microcracked  values.  For 
this  reason,  the  properties  of  micro-cracked  materials  are  expected  to 
exhibit  a strong  dependence  on  temperature  as  well  as  on  thermal  history. 

This  aspect  should  be  taken  into  account  in  assessing  the  thermal  stress 
resistance  of  a micro-cracked  material  intended  for  a specific  engineering 

I 

application. 
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In  summary,  expressions  were  presented  for  the  effect  of  micro- 
cracking on  thermal  conductivity.  Young's  modulus,  strength, Poisson ' s 
ratio  and  the  coefficient  of  thermal  expansion.  The  results  indicate 
that  micro-cracked  materials  should  exhibit  excellent  thermal  stress 
resistance.  Oriented  micro-cracks,  in  principle,  can  be  used  to  consider- 
able advantage.  Also,  micro-cracked  materials  can  be  very  useful  for 
purposes  which  require  high  thermal  stress  resistance  in  combination  with 
high  thermal  insulating  properties. 


Acknowledgment 

The  present  study  was  conducted  as  part  of  a larger  research  pro- 
gram on  the  thermal-mechanical  and  thermal  properties  of  structural 
ceramics  supported  by  the  Office  of  Naval  Research. 


References 


1.  R.  C.  Rossi,  "Thermal-Shock  Resistant  Ceramic  Composites,"  Ceramic 
Bulletin, '48  (7)  736-37  (1969). 

2.  D.  P.  H.  Hasselman,  "Unified  Theory  of  Thermal  Shock  Fracture  Ini- 
tiation and  Crack  Propagation  in  Brittle  Ceramics,"  J.  Am.  Ceram. 

Soc.,  52  (11)  600-04  (1969). 

3.  D.  P.  H.  Hasselman,  "Thermal  Stress  Crack  Stability  and  Propagation 
in  Severe  Thermal  Environments,"  pp.  89-103  in  Materials  Science 
Research,  Vol.  V,  Ceramics  in  Severe  Environments,  Ed.  by  W.  W. 

Kriegel  and  Hayne  Palmour  III,  Plenum  Press,  New  York  (1971). 

4.  E.  A.  Bush  and  F.  A.  Hummel,  "High-Temperature  Mechanical  Properties 
of  Ceramic  Materials:  I,  Magnesium  Dititanate,"  J.  Amer.  Ceram.  Soc., 

41  (6)  189-95  (1958). 

5.  E.  A.  Bush  and  F.  A.  Hummel,  "High-Temperature  Mechanical  Properties 
of  Ceramic  Materials:  II,  Beta-Eucryptite,"  J.  Amer.  Ceram.  Soc.,  42 
(8)  388-91  (1959). 

6.  J.  A.  Kuszyk  and  R.  C.  Bradt,  "Influence  of  Grain  Size  on  Effects  of 
Thermal  Expansion  Anisotropy  in  MgT^O^,"  J.  Amer.  Ceram.  Soc.,  56 
(8)  420-23  (1973). 

7.  H.  J.  Siebeneck,  D.  P.  H.  Hasselman,  J.  J.  Cleveland  and  R.  C.  Bradt, 
"Effect  of  Microcracking  on  the  Thermal  Diffusivity  of  Fe2Ti05", 

J.  Amer.  Ceram.  Soc.,  59  (5-6)  241-44  (1976). 

8.  H.  J.  Siebeneck,  J.  J.  Cleveland,  D.  P.  H.  Hasselman  and  R.  C.  Bradt, 
"Grain  Size  Microcracking  Effects  on  the  Thermal  Diffusivity  of  MgTiO^." 

9.  W.  R.  Manning,  G.  E.  Youngblood  and  D.  P.  H.  Hasselman,  "Effect  of 
Microcracking  on  the  Thermal  Diffusivity  of  Polycrystalline  Aluminum 
Niobate",  J.  Am.  Ceram.  Soc.,  60,  469-70  (1977). 

16 

^ ^ ^ ~ — — •-  - ; 


References  (Cont'd) 


10.  J.  N.  Goodier  and  A.  L.  Florence,  pp.  562-68  in  Proceedings  Xlth 
Int.  Congress  of  Appl.  Mechanics,  Munich  (1964). 

11.  J.  R.  Willis,  "Bounds  and  Self-Consistent  Estimates  for  the  Overall 
Properties  of  Anisotropic  Composites,"  J.  Mech.  Phys.  Solids  25,  185- 
202  (1977). 

12.  D.  P.  H.  Hasselman,  "Effects  of  Cracks  on  Thermal  Conductivity,"  J. 

Comp.  Mat.  12  403-07  (1978). 

13.  J.  P.  Berry,  "Some  Kinetic  Considerations  of  the  Griffith  Criterion 

For  Fracture,"  I and  II,  J.  Mech.  Phys.  Solids,  8 (3),  194-206; 

206-17  (1960). 

14.  J.  B.  Walsh,  "Effect  of  Cracks  on  the  Compressibility  of  Rocks,"  J. 

Geophys.  Res.,  70  (2)  381-99  (1965). 

15.  D.  P.  H.  Hasselman,  "Analysis  of  the  Strain  at  Fracture  of  Brittle 
Solids  with  High  Densities  of  Microcracks,"  J.  Am.  Ceram.  Soc.,  52 
(8),  458-59  (1969). 

16.  R.  L.  Salganik,  "Mechanics  of  Bodies  with  Many  Cracks,"  Izv.  An.  SSSR 
Mekhanika  Tverdogo  Tela,  8 (4)  149-158  (1973). 

17.  R.  J.  O'Connell  and  B.  Budiansky,  "Seismic  Velocities  in  Dry  and  Satu- 
rated Cracked  Solids,"  J.  Geophysical  Research,  79  (35)  5412-26  (1974). 

18.  W.  R.  Delameter , G.  Hermann  and  D.  M.  Barnett,  "Weakening  of  an  Elastic 
Solid  by  a Rectangular  Array  of  Cracks,"  J.  Appl.  Mech.  (43),  74-80 
(1975). 

19.  B.  Budiansky  and  R.  J.  O'Connell,  "Elastic  Module  of  a Cracked  Solid,” 
Int.  J.  Solid  Structures,"  12,  81-97  (1976). 

20.  R.  A.  Sack,  "Extension  of  Griffith's  Theory  of  Rupture  to  Three 
Dimensions,"  Proc.  Phys.  Soc.,  London  729-36  (1946). 


References  (Cont'd) 


21.  B.  Paul,  "Prediction  of  Elastic  Constants  of  Multiphase  Materials," 
Trans.  AIME,  218  (2),  36-41  (1960). 

22.  T.  Yokobori  and  M.  Ichikawa,  "Elastic  Solid  with  an  Infinite  Row  of 
Collinear  Cracks  and  the  Fracture  Criterion,"  J.  Phys.  Soc.,  Japan, 
(19)  2341-42  (1964). 

23.  P.  S.  Turner,  "Thermal-Expansion  Stresses  in  Reinforced  Plastics," 

J.  Res.  Natl.  Bur.  Standards,  37  (4)  239-50  (1946):  RP  1745. 

24.  R.  C.  Rossi,  "Thermal  Expansion  of  BeO-SiC  Composites,"  J.  Amer. 

Ceram.  Soc.,  52  (5)  290-91  (1969). 

25.  D.  P.  H.  Hasselman,  "Thermal  Stress  Resistance  Parameters  for  Brittle 
Refractory  Ceramics:  A Compendium,"  Amer.  Ceram.  Soc.  Bull.,  49  (12), 
1933-37  (1970). 

26.  D.  P.  H.  Hasselman  and  G.  E.  Youngblood,  "Enhanced  Thermal  Stress 
Resistance  of  Structural  Ceramics  with  Thermal  Conductivity  Gradient," 
J.  Amer.  Ceram.  Soc.,  61,  (1)  49-52  (1978). 


Figure  2.  Effect  of  randomly  oriented  micro-cracks  on 
Young's  modulus. 


Figure  3.  Effect  of  randomly  oriented  micro— cracks 
on  Poisson’s  ratio. 
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Figure  4.  Flat  plate  with  rectangular  array  of  cracks  (after 
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footnote  18) . 

•*- 

22 

r 

1 

• 

- -J 

2a/c 


Figure  5.  Effect  of  interacting  cracks  on  Young's 
modulus  of  plate  with  rectangular  array 
of  cracks  shown  in  Fig.  4.  (After  footnote 
18). 
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STRESS  INTENSITY  FACTOR 


Figure  6.  Effect  of  crack  interaction  on  stress  intensity  factor  for  plate 
with  rectangular  arrav  of  cracks  shown  in  Fig.  4.  (After  foot- 
note 18) . 
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Figure  7.  Effect  of  micro-cracking  on  coefficient 
of  thermal  expansion  of  composites  (see 
text  for  data) . 
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Figure  8.  Effect  of  randomly  oriented  micro-cracks 

on  the  thermal  stress  resistance  parameters 
R and  R' . 
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An  analysis  is  presented  of  the  role  of  the  physical  properties 
which  affect  the  thermal  buckling  behavior  of  ceramic  materials  for 
various  geometric  and  thermal  conditions.  For  a column  with  slight 
initial  curvature,  or  a curvature  as  the  result  of  a transverse  tem- 
perature gradient,  expressions  are  derived  for  the  maximum  increase 
in  temperature,  AT^  x or  transverse  temperature  gradient  or  heat  flux 
to  which  the  column  can  be  subjected  such  that  failure  under  the  in- 
fluence of  the  bending  stresses  is  avoided.  From  the  solutions  ob- 
tained a number  of  "thermal  buckling  resistance  parameters"  are  defined 
for  the  comparison  of  the  relative  thermal  buckling  resistance  of 
brittle  materials.  Numerical  examples  indicate  that  thermal  buckling 
of  brittle  materials  of  the  proper  geometry  can  occur  relatively  easily. 
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I. 


Introduction 


Theoretical  analyses  ^ or  the  thermal  stress  resistance  of 
brittle  ceramics  generally  are  based  on  the  assumption  that  failure 
occurs  in  tension  or  shear.  Although  an  extensive  literature^  ^ on 
the  analysis  of  the  thermal  buckling  in  engineering  structures  is 
available,  little  consideration  appears  to  have  been  given  to  the 
possibility  of  thermal  buckling  of  ceramic  materials  used  for  structu- 
ral purposes. 

£ 

In  previous  theoretical  studies,  Fridman  and  Burgermeister  and 
Steup2  presented  solutions  for  the  thermal  buckling  of  an  initially 
straight  column  with  uniform  cross-section  prevented  from  expanding 
against  rigid  constraints.  Under  these  conditions  the  column  remains 
straight  up  to  a critical  temperature  difference,  AT  at  which  the  column 
exhibits  instability  which  (for  column  with  ends  free  to  rotate)  is 
given  by: 


AT  = 7r2I/L2aA 
c 


(1) 


where  I is  the  cross-sectional  moment  of  inertia,  L is  the  length  of 
the  column,  a is  the  coefficient  of  thermal  expansion  and  A is  the  cross- 
sectional  area  of  the  column. 

For  AT  ^ AT^ , the  column  exhibits  post-thermal  buckling  bending 

g 

analyzed  by  Boley  and  Wiener  . For  the  failure  criterion  that  the  maxi- 
mum bending  stresses  in  the  column  should  not  exceed  the  tensile  strength, 


it  was  shown  recently^  that  the  maximum  temperature  difference,  AT 
over  which  a column  with  square  cross-section  can  be  heated,  is: 


max 


a 2 2 2 2 

AT  = AT  + S/L  /it  otE  d 
max  c t 


(2) 
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where  S is  the  tensile  strength.  E is  Young's  modulus  and  d is  the 
smaller  dimension  of  the  cross-srction . 

For  purpose  of  the  analyses  to  be  presented  in  the  present 
paper,  it  should  be  emphasized  that  the  solution  for  AT^  represents 
an  instability  problem.  In  contrast,  the  solution  for  AT  considers 
failure  in  tension  as  the  result  of  the  bending  moment  induced  as  the 
result  of  the  curvature  of  the  column  after  instability  has  occurred. 

The  purpose  of  this  study  is  to  present  additional  solutions 
for  the  thermal  buckling  of  ceramic  materials  for  geometric  and  ther- 
mal conditions  of  interest  to  ceramic  technology.  As  for  the  afore- 
mentioned studies,  the  specific  specimen  geometry  chosen  for  the  present 
analyses  consists  of  a long  slender  column  (rod  or  narrow  plate)  con- 
strained from  expanding  at  the  ends  so  that,  on  a temperature  increase, 
the  column  (etc.)  is  subjected  to  a compressive  load.  For  this  geometry 
and  thermal  loading  the  buckling  behavior  is  analyzed  for,  a)  a column 
which  exhibits  a slight  initial  curvature,  b)  an  initially  straight  column 
with  a curvature  caused  by  a temperature  gradient  in  the  column  per- 
pendicular to  the  column  length  and  c)  an  initially  curved  column  with 
a temperature  gradient  as  for  b,  which  causes  either  an  additional 
positive  or  negative  curvature.  Under  all  conditions,  the  column  will  be 
considered  to  have  a uniform  cross-section  along  its  total  length.  Further- 
more for  simplicity,  the  analysis  will  be  limited  to  deflections  which 
are  sufficiently  small  that  the  curvature  (R  *)  of  the  column  to  a good 
approximation  can  be  expressed  by: 

-1  9 9 9 _9  / 9 9 9 

R = d'y/dx  [1  + (dy/dx)"]  = d v/dx  (3) 

i.e.,  such  that  the  square  of  the  slope  is  small  compared  to  unity  and 
that  the  total  deflections  are  much  less  than  the  column  thickness. 
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The  incidence  of  buckling  of  a column  is  a function  of  the  mechani- 
cal constraints  at  the  column's  ends.  The  present  analyses  will  be  carried 
out  for  columns  which  are  free  to  rotate  at  their  ends,  with  an  end  geometry 
such  that  on  rotation  the  effective  length  of  the  column  between  the  con- 
straints and  therefore  the  mechanical  load,  P = aEATA  remain  invariant.* 

For  different  boundary  conditions,  solutions  can  be  obtained  by  the  general 

. , . . . , 6,7,8,9,10 

mathematical  techniques  given  in  the  literature.  Differences  in 

boundary  conditions  will  affect  only  the  geometric  constants  in  the  solutions 

obtained,  not  the  relative  role  of  the  physical  properties,  easily  verified 

by  dimensional  analysis. 


II.  Analysis 

A.  Uniformly  Heated  Column  with  Slight  Initial  Curvature. 

For  this  column,  the  deflection  (w)  in  the  y-direction  due  to  a uni- 
form initial  curvature  with  radius  R can  be  expressed: 


w = (x  - Lx)/2R 

(12) 

The  condition  of  equilibrium  for  the  column  is: 


Iy  = -aATA{y  + (x  - Lx)/2R} 


(4) 


(5) 


where  y is  the  deflection  of  the  column  in  addition  to  the  deflection 
due  to  the  initial  curvature. 


* For  a rod  of  finite  width  and  square  ends  held  between  perfectly  rigid 

planar  constraints,  the  load  will  increase  for  deflections  of  the  column  up 
to  about  twice  the  column  thickness.  The  assumption  of  constant  load  implies 
that  the  geometric  conditions  of  the  column  ends  and  constraints  are  such  that 
for  small  deflections,  on  rotation  the  net  length  of  the  column  between  the 
constraints  remains  invariant.  A detailed  analysis  of  such  geometric  condi- 
tions is  beyond  the  scope  and  purpose  of  this  paper.  In  general,  for  other 
boundary  conditions  such  as  large  deflections  or  fixed  ends,  the  non-exten- 
sible  bending  of  the  column  needs  to  be  considered.  This  turns  the  problem 
into  a non-linear  one,  which  is  less  convenient  for  the  derivation  of  the 
thermal  buckling  resistance  parameters. 


Eq.  5 has  the  general  solu'  von: 


y = B sin(kx)  + C cos(kx) 


(x2-  1.x) /2R  + (k2R)  1 


(6) 
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where  k = oiATA/l  and  B and  C are  constants. 

From  the  boundary  conditions  y = 0 when  x = 0,  L and  when 
x=L/2,y=0eq.  6 becomes: 


y = -(Rk2)  ^tan(kL/2)  sin  (kx)  - (Rk‘)  ^ cos  (kx) 

-(x2  - Lx)/2R  + (Rk2)-1  (7) 


The  bending  moment  due  to  the  load  P = aEATA  is  a maximum  at  the 

value  y = y (x  = L/2) . From  eq . 7: 
max 

y = -(Rk2r]{sec(kL/2)-l}+L2/8R  (8) 

max 

The  total  deflection  (y^)  of  the  column  at  x = L/2  is  the  sum 
of  ymax  and  w(x  - L/2)  of  eq.  4,  which  yields: 

yt  = -(Rk2)_1{sec(kL/2)-l}  (9) 

For  the  purpose  of  deriving  thermal  buckling  resistance  parameters 
it  is  more  convenient  to  express  eq.  9 by  an  expansion. 

Substitution  of  the  first  three  terms  of  the  general  series 

sec  z = 1 + z 2 / 2 ! + 5z4/4!+  ...  (10) 

where  z = kL/2.  Dividing  eq.  10  by  Rk^  and  retaining  terms  linear  in 
AT  only  yields: 

yt  = -5k2L4/384R  - L^R  (11) 

The  maximum  tensile  stress  (a  ) in  the  column  is  the  sum  of  the  bend- 

max 

ing  stresses  (a,  ) provided  by  the  moment  of  the  total  deflection  and 
b 


thermal  load  P = aEATA  and  the  compressive  stress  oc 
column . 

Using  ab  = Mc/I 


aEAT  in  the 


(12) 


where  M = Y P and  c is  the  distance  from  the  neutral  axis  to  the  outer 

fiber  of  the  column,  yields  for  the  maximum  tensile  stress,  (a  ): 

1 max 


„ - ?»V)2EA2L4C  + aE1T  {aA  _ 

*“  384RI2  8RI 


(13) 


In  eq.  13,  for  slender  columns  (high  L/d  values)  the  major  contri- 
bution to  the  bending  stresses  is  made  by  the  second  term  linear  in  AT, 
followed  by  the  first  term  quadratic  in  AT.  The  compressive  stresses  are 
expected  to  be  minor  only.  The  bending  stresses  linear  in  AT  arise  from 
the  deflection  of  the  column  due  to  the  initial  curvature  which  is  inde- 
pendent of  AT.  The  bending  stresses  quadratic  in  AT  are  due  to  the  de- 
flections which  occur  in  addition  to  the  original  deflection.  These 
additional  deflections  are  a function  of  AT,  which  results  in  a quadratic 
or  higher  order  dependence  on  AT,  depending  on  the  number  of  terms  re- 
tained in  the  expansion  of  eq.  10. 

A solution  of  AT  , the  maximum  temperature  difference  to  which  the 
max 

column  can  be  heated  without  fracture,  can  be  obtained  from  eq.  13,  by 

substitution  of  the  tensile  strength,  S for  a and  solving  the  following 

C IDflX 

quadratic  equation: 


5a2 (AT) 2 EA2LAc  aEAT  ..2 

max  + max^AL  c ^ 


384RI2S, 


8RI 


(14a) 
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In  view  of  its  quadratic  form,  eq . 14 a is  less  convenient  for  the 
derivation  of  thermal  resistance  parameters.  For  this  purpose  two 
specific,  albeit  unusual  cases,  will  be  examined.  If  in  eq.  13  the 
compressive  stresses  are  equal  to  the  bending  stresses  linear  in  AT, 
eq.  13  becomes: 

a = 5a2E(AT)2A2L4c/384RI2  (14b) 

max 

which  for  o * S^.  yields  for  AT  : 

max  t J max 

AT  = (S  /a2E)1/2(384RI2/5A2L4c)1/2  (14c) 

max  t 

If,  in  eq.  13,  the  compressive  stress  term  is  equal  to  the  bending 
stress  quadratic  in  AT: 

a = ctEATAL2c/8RI  (15a) 

max 

which  yields: 

AT  = (S  /aE)(8RI/AL2c)  (15b) 

max  t 

B.  Straight  column  with  transverse  temperature  gradient. 

The  initially  straight  column  to  be  considered  is  subjected  to  a 
uniform  heat  flow  perpendicular  to  the  length  of  the  column  and  paral- 
lel to  the  direction  of  the  shortest  dimension  of  the  cross-section. 

Such  heat  flow  under  steady-state  conditions  and  temperature  indepen- 
dent properties  will  cause  a linear  temperature  distribution  (constant 
temperature  gradient)  through  the  thickness  of  the  column.  In  general, 
a body  of  any  shape  subjected  to  a uniform  temperature  gradient  V = dT/dy 


g 

will  deform  to  exhibit  a uniform  curvature  with  radius  R equal  to  : 


R = (aV) 


-1 


(16) 


where  a is  the  coefficient  of  thermal  expansion.  As  a result,  the 
analysis  of  the  thermal  buckling  of  a straight  column  with  transverse 
temperature  gradient  in  the  y direction  is  identical  to  the  analysis 
of  a column  with  slight  initial  curvature  considered  under  A,  simply 
by  replacing  R in  eqs.  4 through  16  by  (aV)  3 , which  yields: 


a - + aEAT  -1} 

max  384i2  81 


(17) 


Note  that  in  contrast  to  eq.  13,  the  first  and  second  terms  in  eq.  17 
are  proportional  to  the  cube  and  square  of  the  coefficient  of  thermal 
expansion,  resp. 

Again,  similar  to  eq.13,  AT  for  a given  gradient  V,  can  be  ob- 

max 

tained  from  eq.  17  by  substituting  St  for  °max  and  solving  the  quadratic 
equation: 


5a3 (AT) 2 EVA2LAc 
max 

384lV 


ctEAT 


max 


„„  2 
j-qVAL  c 

1 81 


1} 


(18a) 


For  the  purpose  of  obtaining  the  appropriate  thermal  stress  resis- 
tance parameters,  in  analogy  to  eq.  14a,  it  is  convenient  to  derive  AT  ^ 

and  V for  specific  cases, 
max 

If  in  eq.  18a  the  compressive  stress  equals  the  bending  stress  linear 

in  AT,  AT  for  the  straight  column  subjected  to  a transverse  temperature 
max 

gradient  V becomes: 


AT  = (S  /a3E)1/2(384I2/5A2L4Vc)1/2 
max  t 


For  a given  value  of  AT,  the  maximum  temperature  gradient  V 


(18b) 


becomes : 


V = (S  /ct3E){384I2/5(AT)2A2L4c} 
max  t 


(18c) 


Eq.  18c  also  can  be  expressed  in  terms  of  the  maximum  heat  flux, 

Q = KV  , where  K is  the  thermal  conductivity,  which  yields: 
max  max 


2.2„  4 


Q o = (S  K/aJEj(384r/5(AT)^A4L  c 
max  t 


(I8d) 


If  in  eq.  17  the  compressive  stresses  are  equal  to  the  tensile  bending 

stresses  quadratic  in  AT,  AT  can  be  derived  to  be: 

max 


AT  = (S„/ct2E)(8I/AL2Vc) 
max  t 


Similarly,  for  a given  value  of  AT,  V becomes: 

° max 


(19a) 


V = (S„/a2E){8l/AL2(AT)c} 
max  ti 


(19b) 


which  for  the  maximum  transverse  heat  flux,  yields: 


0 = (S  K/cx2E){8I/AL2(AT)c} 

'max  t 


(19c) 


C.  Curved  column  with  transverse  temperature  gradient. 

This  column  with  slight  initial  curvature  of  radius  R is  shown  in 
Fig.  ID,  is  subjected  to  the  identical  thermal  conditions  as  the  column 


considered  under  example  B above.  The  effect  of  the  transverse  heat 
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flow  is  to  add  an  additional  component  Va  to  the  original  curvature  of 
radius  R such  that  the  final  radius  can  be  written: 

Rf  = R(1  + aV) _1  (20) 

In  eq.  20,  it  should  be  noted  the  quantity  7 can  have  a positive  or 
negative  value  depending  on  the  direction  of  heat  flow.  Heat  flow 
directed  towards  the  concave  side  (V  negative)  will  decrease  the  curva- 
ture (increase  R^) , whereas  heat  flow  directed  towards  the  convex  side 
(V  positive)  will  increase  the  curvature.  For  this  column,  the  analysis 
is  identical  to  the  initially  curved  column  or  the  straight  column  with 
a curvature  caused  by  a transverse  heat  flow,  simply  by  replacing  R in 
eqs.  4 through  15  by  R^  as  expressed  by  eq.  20,  with  an  appropriate  change 
in  sign  depending  on  the  direction  of  heat  flow.  Heat  flow,  directed  at 
the  concave  side,  will  cause  a decrease  in  the  maximum  stresses  due  to 
the  reduction  in  the  net  deflection  of  the  column.  Conversely,  heat 
directed  at  the  convex  side  of  the  column,  which  increases  the  net  de- 
flection, will  lead  to  an  increase  in  the  value  of  bending-stresses  and 
a corresponding  reduction  in  the  AT 

max 

An  interesting  condition  occurs  when  the  decrease  in  curvature 
due  to  heat  flow  at  the  concave  side  is  identical  to  the  original 
curvature  of  the  column.  This  occurs  when: 

|V|  = (aR)"1  (21) 

Under  this  condition,  the  column  is  straight  and  the  compressive  load 
caused  by  the  constraints  causes  no  bending  moment  and  corresponding 
bending  stresses  but  places  the  column  in  pure  compression  with  a stress 
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- - - 


of  ctEAT.  Rather  than  a bending  problem,  the  column  is  governed  by 
the  criterion  of  stability  for  a straight  column  as  expressed  by  eq.  1. 
For  values  of  V in  excess  of  the  value  given  to  eq.  22,  the  column 
will  take  on  a net  curvature  opposite  in  direction  to  the  initial  curva- 
ture in  the  absence  of  heat  flow. 

III.  Discussion,  Numerical  Examples  and  Conclusions 

From  eq.  14 , 15,  18  and  19  a number  of  new  thermal  stress  resistance 
parameters  can  be  obtained  appropriate  to  thermal  buckling.  Together 
with  those  from  eqs.  1 and  2 presented  previously^ , these  parameters 
can  be  defined  by: 

Rlb  ' »'1  R2b  * <St2/“E2>  R3b  ' <VaE> 


= (St/u2E)1/2 

R5b  ' (V°2e) 

R,  = (S  K/a2E) 

6b  t 

(22) 

= (St/a3E)1/2 

RSb  * <V“3e) 

R9b  = <StK/a3E) 

It  may  be  noted  that  the  role  of  the  physical  properties  in  these 
parameters  differs  from  that  in  parameters  for  other  thermal  environments 
and  criteria  for  thermal  stress  resistance^  With  those  given  by 

eq.  22  this  brings  the  number  of  non-redundant  thermal  stress  resistance 
parameters  to  a total  of  some  twenty-five.  This  large  number  again  serves 
to  emphasize  the  well-known  conclusion  that  the  relative  thermal  stress 
resistance  of  different  ceramic  materials  can  be  established  only  if  the 
thermal  conditions  are  well  known  quantitatively  and  the  performance 
criteria  and  mode  of  failure  are  pre-determined . 

It  should  be  noted  that  in  the  general  case  of  buckling  failure  as 

described  by  eqs.  13  and  17,  AT  , V or  Q are  described  by  at 

max  max  max 
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least  two  of  the  parameters  given  in  eq.  22,  their  relative  influence 

decided  by  the  geometry  of  the  column. 

Note  that  the  parameter  R, , depends  only  on  the  coefficient  of 

lb 

thermal  expansion  and  in  contrast  to  the  other  parameters  does  not 
depend  on  strength  and  Young's  modulus  as  variables.  This  difference 
is  pointed  out  to  re-emphasize  that  eq.  1 defines  a condition  of  insta- 
bility only  for  a straight  column  and  not  one  of  failure  which  is  deter- 
mined by  the  post-thermal  buckling  behavior  defined  by  eq.  2 and  the 
parameter  R2b'  ^he  seven  remaining  parameters  were  derived  on  the  basis 
of  a definition  of  fracture,  rather  than  instability  and  therefore 
include  strength  and  Young's  modulus  as  variables.  It  may  also  be  noted 

that  with  the  exception  of  R_,  that  the  parameters  of  eq.  22  can  be  ex- 
lb 

pressed  in  fracture-mechanical  terms,  by  means  of  the  appropriate  rela- 
tions between  strength,  critical  stress  intensity  factor  and  the  crack- 
depth.  Furthermore,  it  may  be  noted  that  eqs.  1,  2,  15,  18  and  19  do  not 

depend  on  the  width  of  the  column.  For  columns  which  are  relatively  wide, 

13 

the  stiffness  in  bending  is  greater  than  for  a thin  column  . This  can 

2 

be  accounted  for  by  replacing  E by  E/(l-v  ),  where  v is  Poisson's  ratio 

and  the  appropriate  substitutions  in  the  parameters  of  eq.  22. 

It  is  of  interest  to  illustrate  the  various  expressions  for  the 

thermal  buckling  resistance  by  a number  of  numerical  examples.  For 

an  initially  straight  column  this  was  done  earlier  \ but  will  be 

included  again  for  comparative  purposes. 

An  initially  straight  column  with  square  cross-section  (with  A = bd 

and  I = bd^L2,  where  b = width  and  d = thickness)  can  be  considered  to  be 

—6  "1  5 -2 

composed  of  an  alumina  substrate  material  with  a - 7x10  °C  , E = 4.1x10  MN.m  , 
St  = 4.1xlO^MN.m  ‘,  d = 0.2  cm  and  L = 20  cms.  Eqs.  1 and  2 yields: 
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92°C 


(23) 


AT  = 12°C  and  AT 
c max 

For  the  same  alumina  substrate  with  slight  initial  curvature  with  a 

radius  of  R = 100  cms , eq.  13  for  AT  yields*: 

max 

AT  * 6. 5°C  (24) 

max 

If  the  curvature  of  R = 100  cm  of  the  above  example  were  the  result 

of  transverse  heat  flow  through  an  initially  flat  or  straight  substrate, 

a value  of  AT  - 6.5°C(eq.24)  would  have  been  achieved  at  a temperature 
max  ^ 

gradient  of  approximately  92°C/cm. 

The  analytical  results  obtained  for  the  curved  column  subjected 
to  a transverse  heat  flow,  in  general  are  of  technical  interest  for 
such  structures  as  refractory  arch-roofs  or  walls  of  circular  vertical 
furnaces  such  as  the  blast  or  basic,  oxygen  furnace.  In  general,  it  is 
the  curvature  which  lends  some  degree  of  stability  to  such  structures 
since  it  tends  to  place  the  refractory  under  some  degree  of  compression 
as  the  result  of  either  gravitational  forces  or  thermal  expansion  with 
external  constraints.  If,  however,  such  a structure  is  heated  to  too 
high  internal  temperatures  or  is  heated  too  rapidly  such  that  excessive 
temperature  gradients  are  encountered,  the  net  curvature  of  the  roof  or 
wall  may  well  be  of  opposite  direction.  Under  these  conditions,  such  a 
structure,  unless  well  constrained,  will  have  a tendency  to  buckle  in- 
ward. For  a linear  temperature  gradient,  this  condition  for  eq.  21  occurs 
at  a value  of  critical  temperature  gradient,  given  by: 

Vcr  * (25> 

*This  value  should  not  be  compared  directly  with  ATmax  of  eq.  23  in  view 
of  the  different  boundary  conditions  used  for  the  derivations  of  eq.  2 and 
15.  A solution  for  ATmax  by  means  of  eq.  13  yields  a value  of  ATmax  = 8ttC. 


Taking,  for  example,  a magnesia  brick  wall  with  a - 12x10  ^°C  ^ and 
R = 1000  cms, 

V _ * 80°Ccm_1  (26) 

cr 

Although  the  numerical  examples  presented  above  were  based  on 
arbitararily  selected  dimensions  and  property  values,  the  numerical  re- 
sults nevertheless  indicate  that  thermal  buckling  should  be  considered 
a potential  failure  mechanism  for  brittle  ceramic  materials  and  should 
be  taken  into  account  in  the  design  of  high-temperature  structures. 
Buckling  can  be  minimized  if  the  structures  are  free  tp  expand  or  are 
constrained  in  a manner  that  either  deflections  are  minimized  or  very 
large  deflections  can  occur  at  low  levels  of  tensile  bending  stresses 
(i.e.  very  high  L/d  ratio),  if  practical. 
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THERMAL  STRESS  RESISTANCE  PARAMETERS 
FOR  BRITTLE  MATERIALS  SUBJECTED  TO 
THERMAL  STRESS  FATIGUE 
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W.  ZdaniewskL 
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Brittle  materials  for  high-temperature  applications  usually  are 
highly  susceptible  to  premature  failure  under  conditions  in  which  they 
are  subjected  to  thermal  stress.  The  selection  of  the  optimum  material 
for  a specific  thermal  environment  and  criterion  of  thermal  stress  resis- 
tance can  be  based  on  so-called  "thermal  stress  resistance"  figures-of- 
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merit  or  parameters.  A number  of  such  parameters  are  available  for 
steady-state  or  transient  heat  transfer  (including  inertial  effects) , 
visco-elastic  effects,  thermal  buckling  as  well  as  the  effect  of  catas- 
trophic crack  propagation  on  physical  properties  of  materials  subjected 
to  levels  of  thermal  stress  in  excess  of  the  failure  stress. 

It  is  conceivable  that  materials  undergo  fatigue  under  conditions 
of  thermal  stress  as  the  result  of  sub-critical  crack  growth  at  stress 
levels  below  the  failure  stress.  This  condition  can  arise,  for  instance, 
in  hollow  tubes  conducting  heat  such  as  furnace  tubes  or  those  used 
in  heat  exchangers  or  other  applications.  It  is  the  purpose  of  this 
note  to  derive  on  the  basis  of  fracture-mechanical  principles,  ex- 
pressions for  the  time-to-failure  of  materials  subjected  to  steady- 
state  thermal  stress,  from  which  appropriate  "thermal  fatigue"  resistance 
parameters  can  be  obtained. 

A material  is  considered  which  is  subjected  to  a maximum  value  of 
steady-state  tensile  thermal  stress,  o : 


o = CaEAT/ (1-v) 


where  a is  the  coefficient  of  thermal  expansion,  E is  Young's  modulus, 

v is  Poisson's  ratio  and  C is  a geometric  constant,  which  for  a circular 

hollow-cylinder  with  a temperature  difference  AT  across  the  wall  in  radial 

9 

heat  flow  can  be  obtained  from  the  study  by  Kent. 


I 


Fatigue  of  the  material  occurs  as  the  result  of  sub-critical  growth 
of  a crack  in  the  mode  I crack  displacement  at  a crack  velocity,  V: 

V = AKj11  exp  (-Q/RT)  (2) 

where  A and  n are  numerical  constants,  is  the  mode  I stress  intensity 
factor,  Q is  the  activation  energy  for  crack  growth,  R is  the  gas  con- 
stant and  T is  the  absolute  temperature.  It  will  be  assumed  that  over 
the  total  extent  of  crack-growth  till  the  time-of-f ailure  the  size  of 
the  crack  relative  to  the  structure  size  is  such  that  the  compliance  of 
the  structure  is  not  affected.  Under  this  assumption  the  stress-inten- 
sity factor  Kj  is  related  to  the  stress  and  crack  size  by 

KI  = CTtYal/2  (3) 

where  Y is  a constant  related  to  the  crack  geometry  and  a is  the  crack 
size. 

As  shown  by  Evans10  the  time-to-failure  (tf)  can  be  expressed: 


Cf  = 


/ 


Ic  KjdKj 


2 2 
CTtY  KIi 


(4) 


where  K is  the  critical  stress  intensity  factor  and  K is  the  initial 
Ic  It 

stress  intensity  factor  at  the  imposition  of  the  thermal  stress  at  time, 
t = o. 

Substitution  of  eq.  1 and  3 into  eq.  4 followed  by  integration 
yields  the  time-to-failure  for  the  condition  K <<  K^c  and  n > 2 as: 


t 


f 


2(l-v)2  exp (Q/RT) 

2~~2  (n-2 ) 

(CaEAT)/Y^'A(n-2)IC|  ,'n  ' 


*: 


rVTTT’r - - 

iwiauc^. L£. : . 


(5) 


Eq.  5 also  can  be  expressed  in  terms  of  the  time-to-f ailure  for 
the  situation  that  the  value  of  AT  is  the  result  of  a heat  flux,  q re- 
lated to  AT  by  q = C'KAT,  where  C'  is  another  geometric  constant  and 
K is  the  thermal  conductivity.  This  yields  for  t^: 

2{ (l-v)KC*  }2  exp(Q/RT) 

tf  {CotEq  }2Y2A(n-2)K  ^n_2)  (6 

I i 

Separation  of  the  geometric  and  material  properties  in  eqs.  5 and 
6 yield  the  "thermal  fatigue"  resistance  parameters: 


(1-v)  exp(Q/RT) 


Rf  2 2, 

a“E  (n-2)A 


and  R^  - 


In  principle  under  transient  temperature  conditions  such  as 
those  encountered  in  thermal  cycling,  the  number  of  cycles  to  failure 
can  be  derived  at  least  for  those  conditions  that  the  crack-growth  at 


any  instant  occurs  in  a quasi-static  manner.  With  this  assumption  the 
"thermal  cycle  fatigue"  parameters  should  be  identical  to  those  given 


by  eq.  7. 

As  a last  remark  it  is  of  interest  to  point  out  that  with  the  two 
parameters  derived  in  this  note,  the  number  of  independent  non-redundant 
thermal  stress  resistance  parameters  available  to  the  materials  engineer 
amounts  to  a total  of  twenty-two.  It  appears  almost  superfluous  to  say 
that  the  selection  of  the  optimum  material  for  a given  environment  in- 
volving high  levels  of  thermal  stress  requires  detailed  specification 
of  the  performance  criterion  and  the  thermal  environment  to  which  the 
materials  are  to  be  subjected. 
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ON  THE  THERMAL  FRACTURE  OF  ICE 
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Brittle  materials  are  known  to  be  susceptible  to  catastrophic 
failure  induced  by  thermal  stresses.  This  is  well  known  for  brittle 
materials  for  engineering  applications  involving  high  temperature [ 1 ] . 
Thermal  fracture,  however,  also  can  occur  in  materials  at  moderate  or 
low  temperature,  such  as  ice,  which  can  fail  in  a brittle  manner  even 
at  temperatures  close  to  its  melting  point.  Thermal  fracture  of  ice 
can  be  demonstrated  experimentally  by  the  immersion  of  ice-cubes  in 
fluids  for  industrial  or  house-hold  use.  The  resulting  thermal  frac- 
ture frequently  is  accompanied  by  acoustic  emission  easily  perceived 
by  the  human  ear. 

The  low  thermal  stress  resistance  of  ice  can  be  demonstrated  by 
a numerical  example.  Considered  will  be  a piece  of  ice  of  spherical 
geometry  initially  at  -20°C  suddenly  immersed  in  H2O  at  a temperature 
of  25°C.  For  an  estimate  of  the  heat  transfer  coefficient,  it  will  be 
assumed  that  heat  transfer  to  the  ice  occurs  by  laminar  convection. 

The  heat  transfer  coefficient,  h,  can  be  expressed [2 ] . 

h * £ (0.56)  (N  -N  )**  (1 

D pr  gr 

where  K is  the  thermal  conductivity  of  the  water,  D is  the  diameter  of 

the  sphere,  N is  the  Prandtl  number  and  N is  the  Grasshoff  number, 
pr  gr 


1 


From  tabulations[3]  of  the  properties  of  water  the  quantity  N 
approximately: 


N -N  - 2x104D3AT  cm~3.°C  1 (2) 

pr  8r 

where  AT  is  the  temperature  difference  between  the  ice  surface  and  the 
mean  temperature  of  the  water.  With  D = 2.5  cms,  AT  = 45°C,  = 

cal.  cm.  3°C.  3s  3[4],  the  heat  coefficient  becomes, 

h - 1.98x10  ^cal.cm.  2°C  3.s  3.  (3) 

As  shown  by  Crandall  and  Ging[5]  the  maximum  value  of  tensile  thermal 
stress  (a  ) in  a spherical  body  heated  in  a convection  heat  transfer  en- 
vironment can  be  expressed: 


o = [oEAT/(l-v)] [26/5(6+2))  (4) 

m 

where  AT  is  the  initial  temperature  difference  between  the  ice  and  the 
water,  a is  the  coefficient  of  thermal  expansion,  E is  Young's  modulus, 
v is  Poisson's  ratio  and  8 is  the  Biot  number  defined  by  6 = Rh/K  where 
R is  the  sphere  radius,  h is  the  heat  transfer  coefficient  and  K is  the 
thermal  conductivity  of  the  ice.  Substitution  of  the  following  litera- 
ture data  for  ice,  a = 1.13x10  4oC  ^[4],  E - 103<3N.m  2,  v = 0.33  and 
-3  -1  -1  -1 

K = 5.4x10  cal. cm  .°C  .s  .[6],  with  the  above  value  of  h (eq.  3) 
yields : 


o = 1.82xl07N.m  2 (5) 

m 


This  value  can  be  compared  with  a reported [6]  value  of  the  tensile 

6 “2 

strength,  of  ice  with  o^.  = 1.5x10  N.m  , which  clearly  indicates  that 


for  thermal  conditions  and  dimensions  chosen  for  the  present  calculation, 
the  thermal  fracture  of  ice  is  unavoidable.  A reduction  in  size  of  the 
ice  to  a diameter  of  approximately  0.1  cms  reduces  the  stress  to  a value 
comparable  to  the  fracture  stress.  By  raising  the  initial  temperature 
of  the  ice  to  its  melting  point,  thermal  fracture  on  heating  can  be 
avoided  also,  regardless  of  the  nature  of  magnitude  of  the  heat  transfer 
involved . 
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Expressions  are  derived  for  the  effect  of  penny-shaped  cracks  of 
various  orientations  on  the  thermal  conductivity  of  solid  materials. 

INTRODUCTION 


Tessellated  stresses  in  composites  due  to  thermal  expansion  mis- 
matches in  polycrystalline  materials  which  exhibit  thermal  expansion 
anisotropy  can  cause  excessive  micro-cracking!*^  High  densities  of 
cracks  can  result  from  mechanical  failure  such  as  in  thermal  shock  or 
from  craze  formation.  The  presence  of  such  cracks  can  cause  pro- 
found changes  in  strength,  elastic  behavior,  fracture  toughness  and 

[2  3] 

the  thermal  conductivity  and  diffusivity.  ’ Heavily  micro-cracked 

[4  5] 

materials  have  a potential  for  extreme  thermal  stress  resistance.  ’ 

The  effect  of  cracks  on  elastic  behavior  has  received  theoretical 
attention. ^ Solutions  for  the  effect  of  cracks  on  the  thermal 
conductivity  in  terms  of  crack  density  and  size  do  not,  however,  ap- 
pear to  be  available.  Such  solutions  are  presented  in  this  paper. 


ANALYSIS 

The  present  analysis  for  the  effect  of  cracks  on  thermal  conducti- 
vity is  based  on  the  assumption  that  cracks  can  be  considered  to  be 
pores  of  very  flat  shape.  Specifically,  the  crack  geometry  is  taken 
to  be  equivalent  to  ellipsoids  of  revolution  for  which  the  minor  axis 
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approaches  zero.  Heat  transfer  by  radiation  or  convection  across  the 
cracks  is  assumed  to  be  absent. 

By  treating  the  cracks  as  flat  pores,  the  thermal  conductivity  of 

materials  with  cracks  can  be  obtained  from  literature  solutions  for 

the  conductivity  of  matrices  with  dilute  dispersions  of  ellipsoidal 
rg  q I 

„i  „„  1 » ‘ xn  order  that  the  thermal  conductivity  of  the  material 


particles. 


with  cracks  can  be  treated  as  a continuum  property,  the  size  of  the 

cracks  will  be  considered  to  be  much  smaller  than  the  size  or  volume 

of  the  matrix  material  under  consideration.  Also,  the  cracks  will  be 

assumed  to  be  distributed  uniformly  throughout  the  matrix  material. 

The  cited  literature  equations  on  which  the  present  analysis  is  based 

will  be  given  directly  in  terms  of  the  thermal  conductivity  of  a matrix 

containing  a second  phase  of  thermal  conductivity  equal  to  zero. 

[12] 

From  the  general  solution  of  Fricke,  the  thermal  conductivity 

(K)  of  a matrix  material  containing  dispersions  of  ellipsoidal  shape 
of  zero  thermal  conductivity  can  be  written: ^>9] 


- K0  - g (T^V 

where  K is  the  thermal  conductivity  of  the  matrix  phase,  a,b,  and  c are 
o 

the  axis  of  ellipsoid,  V is  the  volume  fraction  of  the  dispersed  phase 
and:[9’10l 


x,  = l = / _ - u; 

i abc  o (^2)4+32^2)  (s+c2) 

Closed  form  solutions  for  L of  eq.  2 can  be  obtained  for  oblate 

ellipsoids  (a<b=c).  For  the  present  study  the  crack  geometry  will  be 
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assumed  to  be  of  the  form  of  oblate  ellipsoids,  with  a-*o.  Cracks  with 
three  different  orientations  will  be  considered  as  follows: 

A.  Randomly  oriented  cracks. 

Solution  of  eq.  2 for  randomly  oriented  pores  of  oblate  ellipsoidal 
shape  followed  by  substitution  in  eq.  1,  yields  for  the  thermal  conduc- 
tivity: 

K = Ko  3(1-V) [M~  + (3) 


where : 


w ,0  - ^ sin20s  „ 
M = ( — ) cos  0 

sinJ0 


in  which  cos  0 = a/b. 

For  cracks  with  extreme  ellipsoidal  shape,  as  a/b  ■>  o,  0 -*■  tt/2, 
which  yields: 


M = ira/2b 


Substitution  of  eq.  5 for  M in  eq.  3 and  noting  that  M <<  2,  K/M  >> 
2K  and  V <<  1,  upon  rearrangement  eq.  3 becomes: 


K = K (1  + 2bV/3ra) 
o 


Note  that  as  a -+  o,  and  V -+  o for  infinitely  thin  cracks,  the.  second 

terms  in  eq.  6 becomes  indeterminate.  However,  this  indetermacy  can  be 

2 

removed  by  noting  that  for  an  oblate  ellipsoid,  the  volume  equals  Arcab  /3. 
The  volume  fraction  V occupied  by  N cracks  of  equal  size  per  unit  volume 
becomes : 


V = 4nrNab  /3 
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Substitution  of  eq.  7 into  eq.  6 yields: 

K = K (1  + 8Nb3/9)_1  (8) 

o 

which  gives  the  effect  of  cracks  on  thermal  conductivity  in  terms  of 
the  crack  density  and  crack  size,  rather  than  volume  fraction.  Note 
that  b now  corresponds  to  the  radius  of  the  penny-shaped  crack. 

B.  Oriented  cracks  with  crack-plane  perpendicular  to  direction 
of  heat  flow. 

This  case  corresponds  to  oblate  ellipsoidal  pores  with  the  a-axis 
parallel  to  the  heat  flow.  Solution  of  eqs.  1 and  2 yields: 

K = K - VK/M(1-V)  (9) 

o 

which  with  the  aid  of  eq.  5 results  in: 

K = K <l+2bV/-na)_1  (10) 

o 

Eq.  10  with  the  aid  of  eqs.  5 and  7 and  assuming  V <<  1,  yields: 

K = K (l+8Nb3/3)'"1  (11) 

o 

C.  Cracks  with  heat  flow  parallel  to  plane  of  cracks. 

This  case  corresponds  to  ellipsoidal  pores  with  the  major  axis 
(b  or  c)  oriented  parallel  to  the  heat  flow. 

Solution  of  eqs.  1 and  2 gives: 

K = k - — r— 1 (12) 

K o 1-V  2-M  U ' 

Since,  as  a/b  -*■  o,  M <<  2.  Again,  using  the  approximation  V <<  1, 
eq.  12  yields: 

K = K (13) 

o 
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DISCUSSION  AND  CONCLUSIONS 

The  expressions  for  the  effect  of  cracks  on  thermal  conductivity, 
presented  herein  are  expressed  in  terms  of  crack  density  and  crack  size. 
This  should  prove  to  be  of  practical  value  since  the  density  and  crack 
size  probably  are  determined  experimentally  far  more  easily  than  crack 
width  (i.e.  minor  axis)  and  crack  volume.  Also,  a measure  of  crack 
density  and  crack  size  may  be  obtainable  by  measurements  of  the  effect 
of  cracks  on  elastic  properties,  for  which  solutions  are  available. 

For  instance,  the  effect  of  randomly  oriented  cracks  on  bulk  modulus, 

B can  be  expressed: ^ 

B = B [1  + 16(l-v2)Nb3/9(l-2v)]”1  (14) 

o 

Comparison  with  eq.  8 shows  that  the  relative  effect  of  randomly  orient- 
ed cracks  on  bulk  modulus,  apart  from  the  factors  involving  Poisson's 
ratio,  is  "similar"  to  the  effect  of  randomly  oriented  cracks  on  thermal 
conductivity. 

A comparison  of  the  expressions  for  the  effect  of  cracks  of  various 
orientations  shows  that  the  maximum  reduction  in  thermal  conductivity  is 
achieved  when  all  cracks  are  oriented  perpendicular  to  the  direction  of 
heat  flow.  On  the  other  hand,  cracks  oriented  parallel  to  the  direction 
of  heat  flow  have  no  effect  on  thermal  conductivity.  As  expected,  the 
relative  effect  of  randomly  oriented  cracks  on  thermal  conductivity  lies 
between  the  relative  effect  on  thermal  conductivity  of  cracks  oriented 
perpendicular  and  parallel  to  the  heat  flow. 

The  present  analytical  results  can  be  compared  with  the  experimental 
observation  that  excessive  micro-cracking  in  polycrystalline  iron-tita- 
nate,  which  exhibits  a large  degree  of  thermal  expansion  anisotropy. 


-as-- BBii 
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reduces  the  thermal  conductivity  by  a factor  approximately  equal  to 

three  [3].  If  the  grains  are  considered  to  be  cubical  in  shape,  and 

to  contain  at  least  one  crack  each  of  a radius  equal  to  the  grain  size, 

-3  3 

the  crack  density  N = i , such  that  the  product  N£  equals  unity. 

For  randomly  oriented  cracks,  eq.  8 predicts: 

K « 9Ko/17  (15) 

in  good  agreement  with  observation. 

The  analytical  results  presented  herein  are  expected  to  be  ap- 
propriate to  the  areas  of  rock-mechanics,  composites  and  materials  of 
construction  for  high-temperature  applications,  such  as  refractory 
liners  which  often  are  expected  to  perform  satisfactorily  in  a heavily 
cracked  condition,  as  the  result  of  failure  under  severe  thermal  shock. 
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Reaction-sintered  silicon  nitride  (RSSN) , in  view  of  its  excellent 
creep  resistance,  low  density  and  good  thermal  shock  properties,  is  a 
candidate  material  for  high  temperature  turbine  applications.  The  thermo- 
dynamics and  kinetics  of  the  oxidation  of  RSSN  have  received  considerable 
1 2 

attention  ’ . Depending  on  the  temperature  and  specimen  history,  oxidation 

3 4 5 

can  lead  to  an  enhancement  or  decrease  of  strength  ’ ’ . It  is  expected 

that  oxidation  also  should  affect  the  thermal  diffusivity  and  conductivity 

of  RSSN.  The  present  note  reports  experimental  data  for  the  effect  of 

oxidation  of  RSSN  on  its  thermal  diffusivity. 

Specimens  of  RSSN  in  the  form  of  circular  discs  12.5  mm  in  diameter 

-3 

by  2.5  mm  thick  with  a density  of  approximately  2.2  g.cm  were  prepared 
by  nitriding  of  cold-pressed  silicon  compacts  at  a temperature  of  1400°C, 
in  a mixture  of  90  vol.%  Nj  and  10  vol.%  H9,  following  the  firing  schedule 
described  by  Heinrich^.  The  resulting  silicon  nitride  consisted  of  approxi- 
mately 80%  of  the  o-phase,  the  rest  consisting  of  the  B-polymorph.  Oxidation 
was  carried  out  isothermally  in  air  over  a range  of  temperatures  up  to 
1450°C  for  a period  of  50  hours. 

The  thermal  diffusivity  was  measured  by  the  laser-flash  technique^ 

g 

using  equipment  described  in  detail  elsewhere  . 

Figure  1 shows  typical  data  for  the  thermal  diffusivity  as  a function 
of  temperature  for  the  RSSN  oxidized  at  three  different  temperatures. 

Figure  2 shows  the  thermal  diffusivity  at  a number  of  temperatures  as  a 
function  of  oxidation  temperature,  together  with  the  percentage  wt.  gain. 
Figures  3a  and  3b  show  SEM  micrographs  of  fracture  surfaces  of  the  as- 
received  RSSN  and  oxidized  at  1200°C,  which  indicate  that  the  oxidization 
appears  to  have  in  part  closed  the  originally  open  pore  structure.  X-ray 


analysis  shows  that  oxidation  at  1200°C  for  50  hours  resulted  in  an  ap- 
proximately 20%  conversion  of  the  RSSN  into  a mixture  of  crystobalite  and 

amorphous  silica,  while  oxidation  at  1400°C  for  50  hours  produced  smaller 

2 5 

quantities  of  silica,  in  agreement  with  literature  data  ’ . 

Figure  2 indicates  that  the  effect  of  oxidation  on  the  thermal  dif- 

fusivity  of  RSSN  is  a complex  function  of  oxidation  treatment.  Most 

significant  is  the  large  decrease  in  the  thermal  diffusivity  over  a 

rather  narrow  range  of  oxidation  temperature  near  1100°  and  1200°C. 

Also,  the  thermal  diffusivity  of  the  RSSN  oxidized  at  1200°C  is  nearly 

independent  of  temperature,  a behavior  characteristic  of  vitreous  rather 

than  crystalline  materials.  This  general  behavior  can  be  explained  on  the 

1 2 

basis  of  the  known  general  oxidation  behavior  of  RSSN  ’ . At  temperatures 
up  to  about  1000°C  the  oxidation  rate  slowly  increases  with  temperature 
and  occurs  throughout  the  open  porosity  in  the  RSSN.  At  temperatures 
higher  than  1200°C  the  oxidation  product  forms  a protective  layer  on  the 
surface  which  seals  the  pores  and  inhibits  further  oxidation,  with  a 
corresponding  smaller  decrease  in  thermal  diffusivity  than  for  1200°C.  In 
fact,  at  1400°  this  protective  layer,  in  fact,  appears  to  enhance  the  ther- 
mal diffusivity.  The  apparent  decrease  in  thermal  diffusivity  above  1400°C, 
possibly  may  be  due  to  enhanced  oxidation  as  the  result  of  enhanced  thermal- 
ly activated  diffusion  of  oxygen  through  the  protective  oxide  layer.  The 
relative  independence  of  the  thermal  diffusivity  from  temperature  for  the 
RSSN  oxidized  at  1200°C  suggests  that  the  amorphous  silica  phase  may  have 
formed  a continuous  phase  by  preferential  oxidation  at  grain  boundaries. 

9 

The  specific  heat  of  silica  and  silicon  nitride  have  similar  values  . 
Furthermore,  it  should  be  noted  that  the  wt.  gain  of  the  RSSN  on  oxidation 


is  small  (<2%) . For  these  reasons,  most  of  the  observed  changes  in  the 
thermal  diffusivity  due  to  the  oxidation  are  attributable  to  changes  in 
the  thermal  conductivity. 

A decrease  in  thermal  conductivity  is  critical  to  the  thermal  stress 
resistance  of  RSSN.  The  significant  decrease  in  thermal  conductivity, 
especially  for  oxidation  near  1100°  to  1200°C,  may  well  cause  a measurable 
decrease  in  thermal  stress  resistance.  This  effect  will  depend  on  the 
nature  of  thermal  shock  encountered,  but  is  expected  to  be  most  severe 
during  transient  heating  from  low  ambient  temperatures,  where  the  effect 
of  oxidation  on  thermal  conductivity  is  most  pronounced. 

The  present  data  suggest  that  from  the  point  of  view  of  changes  in 
the  thermal  conductivity  and  diffusivity,  oxidation  of  RSSN  should  be 
kept  to  a minimum.  This  may  be  accomplished  by  the  development  of  high  den- 
sity products  by  improved  nitriding,  impregnation  or  closing  of  pores  by 
vapor  deposition.  Alternatively,  on  prior  oxidation  at  1400°C  should  in- 
hibit the  extensive  oxidation  at  lower  temperatures. 
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Effect  of  temperature  on  the  thermal  diffuslvlty  of  reaction- 
sintered  silicon  nitride  subjected  to  oxidation  in  air. 
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?ig.  2.  Thermal  diffusivity  of  reaction  sintered  silicon  nitride  oxidized 
for  50  hrs.  in  air  as  a function  of  oxidation  temperature. 


Fig.  3.  Scanning  electron  fractograph  of  reaction  sintered 


silicon  nitride,  a,  non-oxidized  and  b,  oxidized 
at  1200°C  for  50  hours  (500x) . 
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EFFECT  OF  CRYSTALLIZATION  ON  THE 
THERMAL  DIFFUSIVITY  OF  A CORDIERITE 
GLASS-CERAMIC 

K.  Chyung,  G.  E.  Youngblood  and  D.  P.  H.  Hasselman 

The  magnitude  of  the  thermal  conductivity  and  diffusivity  of 
solids  is  controlled  by  the  specific  heat,  and  transport  properties 
of  the  carriers  responsible  for  the  heat  transfer  such  as  electrons, 
phonons  and  photons1.  Dielectric  glassy  (non-crystalline)  materials 
with  low  values  of  phonon  mean  free  path  at  low  and  moderate  tempera- 
tures at  which  phonon  transport  is  the  primary  mechanism  for  heat 
transport,  are  expected  to  exhibit  much  lower  thermal  conductivity 
and  diffusivity  than  crystalline  materials  with  much  higher  values 
of  mean  free  path.  For  this  reason,  crystallization  of  a glass  to 
form  a glass-ceramic  is  expected  to  cause  a significant  increase  in 
the  thermal  conductivity  and  diffusivity,  as  confirmed  by  experimental 
data  for  the  effect  of  crystallization  on  the  thermal  diffusivity  of  a 
mica-glass-ceramic^. 

In  general,  the  relative  effect  of  partial  crystallization  of 
a glass  on  the  heat  the  heat  transfer  properties  is  expected  to  be 
a function  of  relative  amount  of  the  crystalline  phase,  its  crystal 
structure  and  nature  of  the  atomic  bonding  and  other  properties  rele- 
vant to  the  transport  of  thermal  energy.  The  phonon  mean  free  path 
and  velocity  are  expected  to  be  much  smaller  in  a mica-crystalline 
phase,  with  Vanderwaals  bonding  between  the  basal  planes  than  in  a 
crystalline  phase  with  primarily  ionic  and  co-valent  atomic  bonds. 

The  latter  type  of  crystalline  phase  is  expected  to  cause  a greater 


relative  effect  on  the  thermal  diffusivity  of  a glass-ceramic  than 
a mica  crystalline  phase.  In  the  present  note  this  hypothesis  is 
verified  by  comparing  experimental  data  for  the  thermal  diffusivity 

of  a cordierite  glass  ceramic  with  corresponding  data  for  a mica 

2 

glass-ceramic  reported  previously  . 

A series  of  samples  1"  high  by  1/2"  in  diameter  of  the  original 
glass  of  a commercial  glass-ceramic*  were  crystallized  by  heating 
to  820®C  in  approximately  45  minutes  followed  by  holding  for  2 hours 
for  crystallite  nucleation.  The  samples  were  then  heated  rapidly 
to  and  held  for  8 hours  at  a range  of  temperature  levels  in  order 
to  promote  crystallite  growth.  The  crystallite  morphology  was  stud- 
ied by  replication  transmission  electron  microscopy.  The  nature  of 
the  crystalline  phase  was  measured  by  x-ray  analysis.  The  thermal 
diffusivity  was  measured  from  room  temperature  to  approximately  750°C 

by  the  laser-flash  technique  with  the  same  equipment  and  procedures 
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followed  for  the  mica  glass-ceramic  . 

Figure  1 shows  the  electron  micrographs  for  four  samples  sub- 
jected to  different  heat  treatments.  Table  I lists  the  density  for 
all  samples  and  the  nature  and  amount  of  crystallinity  for  the  samples 
shown  in  Fig.  1.  Figure  2 shows  the  experimental  data  for  the  thermal 
diffusivity.  No  difference  was  noted  between  the  data  obtained  on 
heating  and  cooling  of  two  of  the  samples  for  which  this  effect  was 
investigated . 

Figure  2 indicates  that  near  room  temperature  the  crystallization 
caused  an  increase  in  the  thermal  diffusivity  by  a factor  in  excess 
of  three.  At  the  higher  temperature  the  relative  increase  is  less. 


Corning  Glassworks,  C9606 


Qualitatively,  these  relative  differences  in  thermal  diffusivity  are 
attributable  to  the  combined  effect  of  the  decrease  in  thermal  con- 
ductivity due  to  increased  phonon-scattering  in  the  crystalline  phase, 
with  increasing  temperature  and  the  simultaneous  increase  in  thermal 
conductivity  due  to  photon  transport  in  the  glassy  phase.  It  should 
be  noted  that  the  present  data  for  the  thermal  diffusivity  of  the 

fully  crystallized  glass-ceramic  (1260°C)  lie  somewhat  below  the  data 
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collected  by  Toulouklan  et  al  for  the  commercial  cordierite  glass- 
ceramic.  In  part,  this  difference  may  be  attributable  to  the  dif- 
ference in  the  heat  schedule  used  for  the  crystallization  of  the 
commercial  material  and  the  present  samples,  which  may  have  resulted 
in  differences  in  the  relative  amounts  of  the  various  crystalline 
phases . 

Comparison  of  the  present  data  with  those  reported  for  the  mica 
glass-ceramic  shows  that  the  original  glasses  for  both  studies  have 
about  the  same  value  of  the  thermal  diffusivity.  On  crystallization, 
however,  the  increases  in  the  thermal  diffusivity  obtained  for  the 
cordierite  glass-ceramic  are  far  greater  than  for  the  mica  glass- 
ceramic.  This  observation  is  in  qualitative  agreement  with  the  origi- 
nal hypothesis  that,  for  a glass-ceramic  a crystalline  phase  without 
Vanderwaals  bonding  is  expected  to  be  far  more  effective  in  increas- 
ing the  thermal  diffusivity  of  the  glass-ceramic  than  a crystalline 
phase  in  part  held  together  by  Vanderwaals  bonding.  This  conclusion 
also  is  demonstrated  by  the  present  data  which  show  that  at  crystal- 
lization temperatures  of  1000°C  and  above,  the  increase  in  thermal 
diffusivity  is  governed  primarily  by  the  change  in  the  nature  and 
relative  amounts  of  the  various  crystalline  phases  rather  than  the 
decrease  in  the  amount  in  glassy  phase. 
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Exact  quantitative  estimates  of  this  phenomenon  in  terms  of 
lattice-dynamical  theory  are  beyond  the  scope  of  the  present  study. 

In  general,  however,  the  present  data  and  those  obtained  for  the  mica 
glass-ceramics  offer  strong  support  for  the  suggestion  that  the  ther- 
mal diffusivity  of  glass-ceramics  can  be  tailored  to  specific  values 
by  careful  control  over  composition  and  heat  treatment. 
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TABLE  I.  Heat  Treatment,  Type  and  Amount  of  Crystalline  Phases 


and  Densities  of  Heat-Treated  Cordierite  Class-Ceramics. 


Heat  Treatment  Crystalline  Phases  Density 

(Temp  and  Time)  and  Volume  % gms/cc 

As-Received  Glass  - 2.65 

820°C,  2 hrs . - 2.66 

820°C-2hrs;  905°C-8  hrs.  Mg  - Petalite  (60-80%)  2.76 


6 - Quartz  (10-20%) 

Enstatite  (Trace) 

Class  (10-20%) 

820°C-2  hrs;  1000°C-8  hrs.  6 - Quartz  (85-90%)  2.92 

Spinel  (5-10%) 

Enstatite  (5%) 

Glass  (2-5%) 

820°C-2  hrs;  1100°C-8  hrs.  - 2.67 

820°C-2  hrs;  1150°C-8  hrs.  Cordierite  (70-80%)  2.66 

Rutile  (5-10%) 
a - Cristobalite  (10-15%) 

Glass  (~  2%) 

820°C-2  hrs;  1260°C-8  hrs.  Cordierite  (70-80%)  2.60 

MgTi05  (10-15%) 
a - Cristobalite  (10-15%) 


Fig.  1.  Replication  electron  micrographs  of  partially  crystallized 
cordierite  glass-ceramics.  Nucleation  heat  treatment:  2 hrs  at 
820°C  followed  by  crystal  growth  for  8 hrs  at  a:  905°C;  b:  1000°C; 
c:  1150°C  and  d:  1260°C.  Bar  indicates  one  ym. 


TEMPERATURE  l*C) 


Figure  2.  Effect  of  crystallization  on  the  thermal 
diffusivity  of  a Cordierite  Glass-Ceramic. 
Open  (solid)  symbols  indicate  measurements 
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Particulate  composites  based  on  mica-dispersions  are  excellent 

candidate  materials  for  many  engineering  applications  in  view  of  their 

1 2 

improved  friction  and  wear  behavior  , enhanced  fracture  toughness  and 

3 4 5 

thermal  shock  resistance  * as  well  as  enhanced  machinability  . The 

mica-phase  of  such  composites  can  also  have  a profound  effect  on  physi- 
cal properties  such  as  strength**,  elastic  behavior**,  thermal  expansion 
and  thermal  conductivity^.  The  present  note  reports  data  for  the  effect 
of  a Ba-mica  dispersed  phase  on  the  thermal  diffusivity  of  polycrystalline 
alumina. 

Samples  of  alumina  containing  up  to  50  vol.%  of  a dispersed  phase 
of  Ba-mica,  identical  in  composition  to  those  studied  previously**,  were 
prepared  by  uniaxial  hot-pressing  mixed  powders  of  nominally  BaMg^A^S^- 
O^qF^  mica*  and  y-A^O^  at  1250°C  and  7,000  psi.  The  mica  was  in  the 
form  of  thin  circular  flakes**  approximately  10  pm  thick  by  30  ym  in  dia- 
meter. The  y-A^O^  had  an  average  particle  diameter  of  1.1  ym.  The 
mica  flakes  were  oriented  with  their  long  dimension  and  basal  plane 
perpendicular  to  the  hot-pressing  direction.  The  densities  of  the 

g 

samples  at  room  temperature  are  listed  in  Table  I.  In  agreement  with 
previous  data,  the  densities  with  increasing  mica-content  showed  an 
approximate  linear  negative  deviation  from  the  theoretical  density  cal- 
culated by  the  rule  of  mixtures,  giving  a total  deviation  of  approximate- 
ly 6%  at  50  vol.%  mica.  The  micro-structures  showed  that  their  devia- 
tions in  density  cannot  be  attributed  to  an  increase  in  pore  phase  with 
increasing  mica  content.  Instead  these  density  variations  are  the  result 
of  the  formation  of  a small  amount  of  a glassy  reaction  product  between 

*Synthane  Taylor,  Mybroy  Ceramics  Division,  Ledgewood,  N.J. 

ICerac  Co.,  Menomee  Falls,  Wisconsin. 
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the  mica  and  alumina  phase.  Possibly,  also,  the  lower  densities  may  be 
due  to  micro-cracking  by  cleavage  within  the  mica  due  to  the  internal 
stresses  which  result  on  cooling  as  the  result  of  the  mismatch  in  the 
coefficients  of  thermal  expansion  of  the  mica  and  alumina. 

The  thermal  diffusivity  parallel  to  the  hot-pressing  direction 

from  25°C  to  approximately  800°C  was  measured  by  the  laser-flash  tech- 
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nique  with  equipment  and  procedures  described  elsewhere  . Fig.  1 
shows  the  experimental  data  for  all  compositions  with  the  exception  of 
the  sample  with  20  vol.%  mica,  which  exhibited  an  apparently  anomalous 
behavior.  In  order  to  provide  a basis  of  comparison  of  the  relative 
temperature  dependence.  Fig.  1 includes  experimental  data  for  a pure 
dense  polycrystalline  alumina*.  Fig.  2 shows  the  compositional  de- 
pendence for  various  temperatures  obtained  from  the  curves  drawn 
through  the  data  of  Fig.  1.  Figure  2 also  includes  the  data  for  the 
sample  containing  20  vol.%  mica. 

The  anomaly  in  the  present  data  for  the  thermal  diffusivity  at 
20  vol.%  mica  was  also  found  for  the  thermal  conductivity^ , with  the 
exception  that  in  the  latter  study,  the  anomaly  occurred  at  15  vol.% 
mica.  This  anomaly  is  attributable  to  the  type  and  amount  of  reaction 
products  formed  between  the  mica  and  the  alumina  during  hot-pressing 
Preparation  of  full-density  composites  is  accomplished  within  a narrow 
temperature  range  below  which  complete  densif ication  does  not  occur 
and  above  which  undesirable  reaction  products  are  formed.  This  situation 
is  particularly  critical  near  15-20%  Ba-mica.  Below  this  amount  an 
excess  temperature  results  in  the  formation  of  a glass  and  a spinel, 

*SN60,  Feldmdhle,  Plochingen,  West  Germany 


the  composite  containing  20%  mica,  smaller  than  for  the  other  samples 
as  expected  from  initial  amount  of  mica  present.  For  the  present 
data  the  anomalous  value  for  the  thermal  diffusivitv  most  likely  is 


because  the  reaction  products  either  have  a higher  thermal  conducti- 


vity or  a lower  specific  heat  than  the  corresponding  value  for  the 
mica  or  the  alumina. 

The  relative  temperature  dependence  of  the  thermal  diffusivity  of 

the  mica-alumina  composites  is  similar  to  that  for  the  pure  alumina. 

In  fact,  a plot  of  a ^ vs  T shows  a straight  line  behavior  with  an 

intercept  near  zero  at  T = 0°K  for  all  compositions.  It  may  be  noted 

again  that  the  present  data  are  for  the  thermal  diffusivity  with  heat 

flow  through  the  mica  occurring  perpendicular  to  the  weakly  bonded  basal 

planes.  As  expected  from  experimental  data  for  other  micas  \ boron 
12  13  14  15 

nitride  ’ and  pyrolytic  graphite  ’ , heat  flow  perpendicular  to 

the  basal  plane  is  expected  to  result  in  a thermal  conductivity  at 
least  an  order  of  magnitude  less  than  that  of  polycrystalline  aluminum 
oxide  and  also  to  result  in  a temperature  dependence  less  than  T , ex- 
pected for  phonon  conduction  over  the  range  of  temperature  of  the  pre- 
sent data.  The  reason  that  the  temperature  dependence  of  the  thermal 
conductivity/dif fusivity  of  the  mica  has  little  or  no  influence  on  the 
temperature  dependence  of  the  Ba-mica/alumina  composites  is  that  the  heat 
is  conducted  primarily  by  the  continuous  alumina  phase.  If,  however, 
the  mica  had  been  the  continuous  phase,  its  temperature  dependence  would 
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have  had  a significant  effect  on  the  temperature  dependence  of  the 
thermal  conductivity  and  diffusivity  of  the  composite,  especially  at 
higher  mica  volume  fractions.  This  latter  conclusion  is  easily  veri- 
fied by  expressions  for  the  thermal  conductivity  of  composites  or  by 
analogy  by  equivalent  parallel  or  series  electrical  circuits.  This 
suggests,  in  general,  that  the  magnitude  of  thermal  conductivity  and 
its  temperature  dependence  in  composites  can  be  controlled  by  the  pro- 
per choice  of  components  and  their  distribution  within  the  composite. 
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TABLE  I.  Composition  and  Density  of  Ba-Mica/Al^O^  Composites. 


Ba-Mica  Content  (Vol.%) 


Density  (gms/cc) 


3.935 


3.845 


3.750 


3.695 


3.625 


3.540 


Fig.  2.  Effect  of  Composition  on  the  Thermal  Diffusi- 
vity  of  Ba-Mica-Alumina  Composites. 
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Hot-pressed  glass-crystal  mixtures  are  used  as  model  materials  for 

studies  of  the  properties  of  brittle  matrix  composites.  Properties  such 

1 2 3 2 4 3 4 

as  elastic  behavior  ’ ’ , strength  ’ and  fracture  toughness  ’ have  been 

studied  extensively  for  glass  matrices  with  an  alumina  dispersed  phase. 
However,  such  composites  in  their  own  right  are  also  candidate  materials 
for  engineering  applications  involving  thermal  stress.  This  note  presents 
data  for  the  thermal  conductivity  and  thermal  diffusivity  of  glass-alumina 
mixtures,  which  permits  making  an  assessment  of  the  effect  of  the  cry- 
stalline phase  on  thermal  stress  resistance. 

The  specific  glass-alumina  specimens  in  this  study  were  prepared 

3 

in  a program  addressing  the  fracture  toughness  of  a borosilicate  glass 

with  alumina  dispersions.  The  glass  consisted  of  70  mole  % Si02  plus 

and  Na20  in  a molar  ratio  of  0.67.  The  alumina  dispersions  were 

spherical  with  a diameter  of  25  ± 7 pm.  Details  for  the  specimen  pre- 

3 

paration  and  microstructure  were  presented  earlier  . This  particular 
glass-alumina  system  was  selected  because  of  the  close  match  between 
the  coefficients  of  thermal  expansion  of  the  two  phases,  in  order  to 
minimize  or  eliminate  the  formation  of  micro-cracks  due  to  internal 
stresses;  such  cracks  could  have  significant  effect  on  the  thermal  con- 
ductivity^ . 

The  thermal  diffusivity  of  the  composites  was  measured  over  the 
temperature  range  from  room  temperature  to  about  600°C  by  the  laser- 
flash  diffusivity  technique^  using  equipment  described  in  detail  else- 
where^, with  the  transient  temperatures  of  the  specimens  during  the 
test  being  monitored  by  IR-detectors.  The  thermal  conductivity  was 
calculated  from  the  thermal  diffusivity  using  values  for  the  density 
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of  the  alumina  and  the  glass  of  3.987  and  2.454  g.cc.  ^ resp.  and 

g 

literature  data  for  the  specific  heat  of  alumina  and  a borosilicate 

9 

glass  of  composition  similar  to  the  glass  of  the  present  study  . 

Figure  1 shows  the  experimental  data  for  the  thermal  diffusivity 
as  a function  of  temperature  for  a number  of  compositions  for  which  the 
volume  fractions  were  determined  from  the  composite  densities  after 
hot-pressing  . For  some  of  these  compositions,  fig.  2 shows  the  cal- 
culated values  of  thermal  conductivity.  The  relative  effect  of  temperature 
and  composition  on  the  thermal  diffusivity  and  conductivity  of  these 
composites  can  be  explained  on  the  basis  of  the  general  heat  conduction 
processes  in  dielectric  materials.  Due  to  the  relative  differences  in 
elastic  properties  and  degree  of  crystallinity,  the  thermal  conductivity 
and  diffusivity  of  the  glass  is  expected  to  exhibit  much  lower  values 
and  a lower  temperature  dependence  than  the  corresponding  values  for  the 
alumina.  For  these  reasons,  addition  of  alumina  particles  to  the  glassy 
phase  is  expected  to  increase  the  magnitude  of  the  thermal  diffusivity 
with  a corresponding  increase  in  its  temperature  dependence,  in  agreement 
with  observations.  The  positive  temperature  dependence  of  the  thermal 
diffusivity  of  the  glass  at  the  higher  temperatures  most  likely  is  attri- 
butable to  a contribution  of  internal  radiation  to  the  heat  transfer  pro- 
cess. 

The  values  of  the  thermal  conductivity  show  a much  lower  relative 
temperature  dependence  than  the  thermal  diffusivity,  because  the  positive 
temperature  dependence  of  the  specific  heat  compensates  for  the  negative 
temperature  dependence  of  the  thermal  diffusivity.  It  may  be  noted  that 
at  the  highest  volume  fraction  the  variation  in  thermal  conductivity  with 
volume  content  alumina  lies  above  the  values  calculated  by  the  Rayleigh- 
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Maxwell  theory^’ which  is  appropriate  only  for  dilute  concentrations 
of  spherical  inclusions. 

An  assessment  of  the  anticipated  effect  of  the  alumina  phase  on 

thermal  stress  resistance  can  be  made  on  the  basis  of  its  influence  on 

thermal  stress  resistance  parameters.  As  a measure  of  the  resistance 

of  a brittle  material  to  the  initiation  of  thermal  fracture  two  well- 
12 

known  parameters  are:  R = St(l-v)/aE  and  R'  = St (l-v)K/aE , in  which 
St  is  the  tensile  strength,  v is  Poisson's  ratio,  a is  the  coefficient 
of  thermal  expansion,  E is  Young's  modulus  and  K is  the  thermal  conduc- 
tivity. The  resistance  of  crack  propagation  and  corresponding  loss  in 
load-bearing  ability  after  thermal  fracture  initiation  described  is  given 
by  the  parameters^,  R = (G/a“E)  2 and  R^t  = (GK^/a^E)^  for  stable 
crack  propagation  and  R""  = GE/St  for  unstable  catastrophic  crack  pro- 
pagation. In  these  latter  parameters  G is  the  fracture  energy,  related 

k 

to  the  critical  stress  intensity  factor,  KT  = (2GE)  . Figure  3 shows 

i c 

the  relative  variation  of  these  parameters  with  volume  fraction  alumina 

3 

as  calculated  from  the  thermal  conductivity  and  mechanical  property  and 
14 

strength  data  for  these  materials.  The  location  of  the  curves  for  the 
parameters  Rgt,  R^t  and  R""  must  be  considered  only  approximate  in  view 
of  the  lack  of  sufficient  data  points  to  plot  the  curves  accurately.  Never- 
theless, sufficient  data  were  available  to  indicate  that  these  parameters 
vary  significantly  with  alumina  content. 

As  judged  by  the  variation  of  the  parameters  R and  R' , the  addition 
of  the  alumina  phase  has  only  a minor  influence  on  the  thermal  conditions 
required  to  initiate  fracture.  In  fact,  the  variation  in  R indicates 
that  the  alumina  phase  actually  decreases  thermal  stress  resistance. 
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The  parameter  R'  increases  primarily  because  the  increase  in  K exceeds 
the  increase  in  E.  As  indicated  by  experimental  data,  the  values  of  R 
and  R'  could  be  increased  by  increasing  the  composite  strength  by  de- 
creasing the  alumina  particle  size. 

The  major  increases  in  the  parameters  R and  R^t  and  especially 

R"",  indicate  that  the  alumina  crystalline  phase  has  a large  effect  on 

the  nature  of  crack  propagation  following  thermal  fracture.  Specifi- 

cally,  under  conditions  of  both  stable  and  unstable  fracture  the  extent 

of  crack  propagation  should  be  reduced  appreciably  accompanied  by  a major 

increase  in  the  retention  of  load-bearing  ability  or  any  other  property 

affected  by  cracks.  The  general  shape  of  the  curves  in  fig.  3 for 

R , R'  and  R""  indicates  that  optimum  resistance  to  thermal  crack 
st  st 

propagation  can  be  obtained  at  approximately  20  vol.%  alumina. 

These  results  suggest  that  by  incorporating  a high  conductivity 
crystalline  phase  in  the  glassy  material,  the  resistance  to  the  initia- 
tion and  the  catastrophic  nature  of  thermal  fracture  of  a glass  can  be 
improved  significantly.  Well  controlled  thermal  stress  experiments  to 
confirm  this  conclusion  should  be  of  interest. 
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ABSTRACT 


A review  is  presented  of  the  selection  rules  for  brittle  structural 
ceramics  subjected  to  severe  thermal  stress.  A total  of  twenty-two  figures- 
of -merit  are  presented  for  a total  of  twenty-nine  different  heat  transfer 
environments  and/or  criteria  of  thermal  stress  resistance,  involving  steady- 
state  or  transient  heat  transfer,  internal  heat  generation,  thermal  buckling 
as  well  as  thermal  fatigue  and  unstable  :rack  propagation. 


1. 


Introduction 


Materials  of  construction  for  high-temperature  structures  and/or 
components,  in  spite  of  such  properties  as  high  melting  points,  low 
vapor  pressure,  high  creep  resistance  and  other  favorable  characteristics, 
tend  to  be  highly  brittle.  As  a direct  result  of  this  brittleness,  high- 
temperature  materials  are  highly  susceptible  to  catastrophic  failure  under 
thermal  conditions  which  involve  appreciable  levels  of  thermal  stress.  In 
order  to  avoid  thermal  stress  failure  it  is  essential  that  the  structure 
or  component  be  designed  to  minimize  the  conditions  which  lead  to  the 
generation  of  thermal  stresses,  and  that  the  material  of  construction  be 
selected  with  a combination  of  properties  such  that  the  magnitude  of  ther- 
mal stress  is  minimized.  In  order  to  achieve  this  latter  objective  it  is 
critical  that  the  role  of  the  pertinent  material  properties  which  govern 
thermal  stress  fracture  be  well  understood.  For  this  purpose,  a large 
number  of  analytical  expressions  have  been  presented  in  the  literature 
devoted  to  materials  science  and  engineering  for  the  thermal  stress  re- 
sistance of  brittle  materials  for  a wide  range  of  thermal  environments  and 
criteria  of  thermal  stress  resistance  in  terms  of  the  appropriate  material 
properties  and  geometric  variables.  From  these  expressions  figures-of- 


merit  or  so-called  thermal  stress  resistance  parameters  can  be  derived 


which  form  the  basis  for  the  selection  of  the  optimum  material  for  a parti- 
cular thermal  environment  and  criterion  of  thermal  stress  resistance. 

These  thermal  environments  include  steady-state  heat  flow,  transient  heat 
flow  involving  convective  or  radiative  heat  transfer,  thermal  fatigue, 
thermal  stability  as  well  as  thermal  stress  resistance  controlled  by  stable 
or  unstable  catastrophic  crack  propagation. 


The  present  writer  some  nine  years  ago  presented  a compendium  of  the 
approximately  ten  or  twelve  thermal  stress  resistance  parameters  available 
at  that  time.  Since  then,  the  total  number  of  non-redundant  equations  for 
thermal  stress  resistance  and  number  of  thermal  stress  resistance  parameters 
has  grown  appreciably.  In  fact,  the  total  number  of  parameters  has  grown 
to  some  twenty-two.  This  indicates  that  the  selection  or  development  of 
material  with  optimum  thermal  stress  resistance  for  a given  thermal  environ- 
ment is  a highly  complex  task.  In  view  of  this  complexity  and  the  gain  in 
understanding  achieved  during  the  last  decade,  it  appears  justified  to 
present  a review  of  the  thermal  stress  f igures-of-nerit  in  order  to  present 
the  reader  with  the  latest  information.  The  presentation  of  such  a review 
represents  the  primary  objective  of  this  paper,  which  is  organized  as 
follows : 

For  a given  thermal  environment  and  criterion  of  thermal  stress 
resistance,  a literature  equation  is  presented  for  a single  geometry. 

On  dimensional  grounds,  for  the  same  thermal  environment  and  criterion  of 
thermal  stress  resistance,  with  exception  of  the  geometric  constant,  the 
equation  presented  also  is  appropriate  for  other  geometries.  For  this 
reason  a general  thermal-stress-resistance  f igure-of-merit  can  be  derived 
appropriate  to  the  specific  thermal  environment  and  performance  criterion 
independent  of  geometric  considerations.  The  appropriate  f igure-of-merit 
is  given  immediately  following  the  equation  given.  Depending  on  the  con- 
dition of  external  constraints  or  geometry,  Poisson's  ratio  can  enter  in 
various  terms  such  as  (1-v) , (l-2v) , (1+v)  etc.,  or  may  even  be  absent. 

For  purposes  of  brevity,  these  terms  involving  Poisson's  ratio  are  con- 
sidered part  of  the  geometric  constant,  which  will  keep  the  total  number 
of  expressions  to  be  presented  to  a manageable  number.  Nevertheless, 


appropriate  literature  references  involving  terms  for  Poisson's  ratio 

will  be  given  whenever  available. 

Again  for  purposes  of  brevity,  the  expression  to  be  presented  wil] 

be  given  with  an  adequate  but  brief  description  of  the  thermal  environment 

or  other  pertinent  conditions.  For  details  of  the  mathematics,  assumptions, 

boundary  conditions,  etc.,  the  reader  is  referred  to  the  original  papers. 

If  an  expression  or  f igure-of-merit  has  not  appeared  in  the  literature 

previously,  it  will  be  so  stated?  Whenever  possible,  the  symbols  used  in 

the  equations  are  those  recommended  by  the  American  Ceramic  Society  and 
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are  listed  separately. 


Equations  and  Figures-of-Merit 


A.  Steady-State  Heat  Flow 

a.  Circular  hollow  cylinder;  radial  heat  flow,  AT  across 

max 

....  3,4 
wall  : 


AT  = S(l-v)C/otE 
max  t 


(la) 


Figure-of-merit : 


R = St(l-v)/«E  (lb) 


b.  In  terms  of  maximum  heat  flow,  0 , eq.  1 becomes: 

max 


Q_  = S. (l-v)KC '/oE 

max  t 

(2a) 

Figure-of-merit : 

R’  = St(l-v)K/aE 

(2b) 

the  equation  presented 

* The  use  of  the  term  "can  be  derived"  implies  that 
can  be  derived  readily  from  the  reference  cited. 

L 

3 

. , 

R'  = S.  (l-v)K/an  (4b) 

cr  t 

e.  Linear  heat  flow  around  spherical  cavity.  Maximum  temperature 
gradient^  can  be  derived  to  be: 

V = 2S. (l-v)/aEb  (5a) 

max  t 

Figure-of -merit : 


R = St(l-v)/aE 

f.  Eq.  5a  in  terms  of  maximum  heat  flow  becomes: 

Q - 2S. (l-v)K/aEb 
max  t 

Figure-of -merit : 


R'  = St(l-v)K/aE 


From  other  solutions  f igures-of-merit  can  be  obtained 
for  inclusions  other  than  pores. 

B.  Internal  Heat  Generation 

g.  For  solid  circular  rod  with  uniform  internal  heat  generation 

g 

(H) , the  maximum  rate  of  heat  generation  is  : 


H = 8S„(l-v)K/c<Eb 
max  t 


Figure-of-merit : 


R'  = St(l-v)K/aE 


C.  Transient  Heat  Transfer 


h.  Flat  plate;  Linear  rate  of  increase  of  surface  temperature. 

9 

Maximum  rate  of  increase  : 


T = 3S„(l-v)a/«Eb 
max  t 


Figure-of-Merit : 


R"  = St(l-v)a/aE 


i.  Convective  (Newtonian)  heat  transfer.  Spherical  body  subjected 

Arr10 

to  heating  and  cooling  over  temperature  range  AT 


2.5  (1-v) 


(1  + 2/8)  8 = bh/K 


For  severe  heat  transfer,  8 -*■  00 


AT  =2.5  S (l-v)/aE 
max  t 


(10a) 


Figure-of-Merit : 


R = St(l-v)/aE 


(10b) 


I 


I 


For  mild  heat  transfer,  3 <<  1 


AT  = 5S  (l-v)K/aEbh 
max  t 


(11a) 


Figure-of-merit : 


R'  = St(l-v)K/aE  (lib) 

j.  Radiation  heat  transfer.  Spherical  body  at  low  initial 
temperature  and  wavelength  independent  emissivity.  In- 
stantaneously uniformly  subjected  to  black-body  radiation. 
Maximum  radiation  temperature^: 


T = (5S. (l-v)K/pbaEe } 
max  t 


(12a) 


Figure-of-merit : 


Rrad  = {St(l-v)K/aEe}4 


(12b) 


k.  Radiation  heat  transfer.  If  body  transparent  to  radiation  at 

12 

wavelength,  X < X0  eq.  12a  becomes  : 


T = {5S. (l-v)K/pb(l-F.  )aEe}4 
max  t Xo 


(13a) 


Figure-of-merit : 


R.  = {S„(l-v)K/(l-F,  )aEe) * 

trans  t Xo 


(13b) 


Transient  heat  transfer  with  inertial  effects.  Rod  of  length 
L,  freely  suspended.  Rod  is  very  rapidly  exposed  to  uniform 
temperature  increase  AT  in  time  period,  tD.  Compressive  stresses 
will  occur  because  of  finite  time  required  for  expansion  to 
take  place. ^ 

The  maximum  temperature  increase  of  the  rod,  without  compressive 
failure  can  be  derived  to  be: 

6 


,u.._ 


M 


AT  = S /«E  for  t It 
max  c o m 


(14a) 


Figure-of -merit : 


R.  = S /aE 
in  c 


(14b) 


AT  = S c t / aEL  for  t t t 

max  co  o m 


(15a) 


Figure-of-merit : 


R'  = S c/aE 
in  c 


(15b) 


in.  Transient  heat  transfer.  Time-to-failure.  Sphere.  Surface 

temperature  raised  instantaneously  by  temperature  difference, 

14 

AT.  Time  to  maximum  stress  : 


t*  = 0.0574  b /a 


06a) 


Figure-of-merit : 


Rf  = a 


(16b) 


D.  Thermal  Buckling 

n.  Straight  column  prevented  from  expanding  against  rigid 

constraint.  Heated  by  temperature  difference  AT.  Critical 
temperature  difference  required  for  column  instability^: 


AT  = n2I/L2Aa 
c 

Figure-of-merit : 


(17a) 


Rlb  = ° 


(17b) 


o.  Post-thermal  buckling  of  straight  column.  Originally  straight 

column  prevented  from  expanding  against  rigid  constraints,  heated 
over  temperature  difference,  AT>ATC  to  be  subjected  to  bending 
in  post-thermal  buckling. 


7 


Maximum  temperature  difference  required  for  failure  of 
square  column  as  the  result  of  tensile  bending  stresses3**: 


22  2 22222 
AT  * it  <T/12L  a + S/L  /t  aE  d 
max  t 


(18a) 


Figure-of-Merit : 


R2b  * St2/aE2 


(18b) 


p.  Column  with  slight  initial  curvature  (radius  R)  prevented 
from  expanding  against  rigid  constraints.  Heated  by  tem- 
perature difference  AT.  Critical  temperature  difference 
required  for  column  fracture  due  to  bending  moment3  : 

AT  = (S./cx2E)^(16  Rd3/15L4)^ 
max  t 


(19a) 


Figure-of-merit : 


R3b  = (St/a  E) 2 


(19b) 


q.  Column  initially  straight  subjected  to  transverse  temperature 
gradient  V and  prevented  from  expanding  against  rigid  con- 
straints. Heated  by  temperature  difference  AT.  Critical 
temperature  difference  for  fracture  of  column  due  to  bending 
moment  (for  given  gradient  V)^: 


VT  = (S. /a  E) z(16d  /15VL  ) 
max  t 


(20a) 


Figure-of-merit : 


R4b  = (St/c3E)% 


(20b) 


r.  Initially  straight  column  subjected  to  transverse  temperature 


gradient  7 and  prevented  from  expanding  against  rigid  con- 
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straints.  Heated  by  temperature  difference  AT.  Critical 
temperature  gradient  at  fixed  AT,  required  for  column  frac- 


ture 

V = (S _/«3E)  {16d3/15(AT)2L4}  (21a) 

max  t 

Figure-of-merit : 

R5b  = (St/a3E)  (21b) 

Eq.  21a  in  terms  of  maximum  heat  flow  becomes3^: 

Q = (S  K/c*3E)  U6d3/15(AT)2L4}  (22a) 

max  t 

Figure-of-merit : 

R,.  = (S  K/a3E)  (22b) 

6b  t 


E,  Crack-Propagation 

s.  Fatigue.  Subcritical  crack  growth.  Constant  thermal  stress. 

Material  subjected  at  t = o to  a steady-state  thermal  stress,  o 

CaEAT/ (1-v) , with  a <<  St>  the  fracture  stress.  The  material 

undergoes  subcritical  crack  growth  with  a crack  velocity,  V = 

AX TexP (-Q/RT) . For  K_  = Ya  faT  and  K_  (t=o)  <<  XT  , for 

I la  1 1 Ic 

the  condition  that  over  the  total  time  period,  the  total  ex- 
tent of  crack  growth  does  not  affect  the  compliance,  the  time 
18 

to  failure  (t^)  is  : 

2(l-v)2  exp  (O/RT)  (23a) 

tf  C(aEAT)2Y2A(n-2)KTi(n'2) 


9 


I 


Figure-of -merit : 


(1-v)  exp  (Q/RT) 
a2E2(n-2)A 


(23b) 


If  in  analogy  to  eq.  2a,  AT  in  eq.  23a  is  related  to  a 

18 

heat  flux  Q,  with  AT  * C'Q/K,  the  time  to  failure  becomes  : 


J2(l-v)2K2  exp  (Q/RT) 
e(aEATC'Q)2Y2A(n-2)K 


(n-2) 


(24a) 


Figure-of-merlt : 


(1-v)2K2  exp  (Q/RT) 
a2E2(n-2)A 


(24b) 


t.  Crack  stability.  Onset  of  catastrophic  crack  propagation- 

Flat  plate  with  crack.  Heat  flow  perpendicular  to  crack. 

19 

Maximum  temperature  gradient  to  avoid  crack  propagation  : 


V = 8(G/2irot2E  i?)'2 
max 


(25a) 


Figure-of-merit : 


R = (G/a  E) 
st 


(25b) 


Eq.  25a  in  terms  of  maximum  heat  flow  becomes: 


Q_  = 8(GK2/2ira2EiV* 

inax 


(26a) 


Figure-of-merit : 


R'  = (GK2/a2E)^ 

St 


(26b) 


u.  Flat  plate  with  N cracks  per  unit  area.  Plate  cooled  by  tem- 
perature difference  AT  and  prevented  from  shrinking  by  uni- 
axial rigid  constraints.  Critical  temperature  difference 


I 


of  cooling 


AT  = (2G/nJla2E)f(l+2TrNi2) 
max 


(27a) 


Figure-of -merit : 


Rst  = (G/a  E) 2 


(27b) 


For  three-dimensionally  constrained  body  with  penny- 
21 

shaped  cracks  AT  yields  same  f igure-of-merit  as  eq.  27b. 
max 

F.  Catastrophic  crack  propagation  and  resulting  change  in  strength  behavior: 

v.  Flat  plate  undergoing  catastrophic  (unstable)  fracture  by 

simultaneous  propagation  of  N cracks  per  unit  area.  The 

20 

final  crack  length  at  crack  arrest  : 


Jtf  = S' /4NGE 


(28a) 


Minimum  extent  of  crack  propagation  requires  high  values  of 


the  f igure-of-merit : 


R""  = GE/S 


(28b) 


The  sane  f igure-of-merit  is  obtained  for  penny-shaped  crack 

21 

undergoing  crack  propagation  in  three-dimensional  body  . In 

2 

order  to  avoid  redundancy,  the  f igure-of-merit , R' ' ' = E/St  , 
proposed  previously  is  not  included  in  the  present  review. 

w.  Strength  retained  (S^)  by  plate  which  has  undergone  crack 

propagation  with  final  crack  length,  ^ can  be  derived  to  be: 


c - £!_  r8V* 
f St  ' 


Figure-of-merit : 


(29a) 


Rf  = GE/St 


11 


(29b) 


Z3CT 


For  other  geometries  eq.  29a  can  take  on  somewhat  different 

, 22 
form 

3.  Discussion  and  Conclusions 

The  expressions  and  f igures-of-merit  presented  above  indicate  that 

the  selection  of  a material  with  optimum  thermal  stress  resistance  can 

be  a complex  task.  A total  of  twenty-nine  independent  equations  are 

given  for  the  thermal  stress  resistance  of  brittle  materials  for  a number 

of  different  thermal  environments,  failure  mechanisms  and  measures 

of  thermal  stress  resistance.  These  twenty-nine  equations  lead  to  a total 

of  twenty-two  different  thermal  stress  figures-of-merit  summarized  in 

Table  I.  It  should  be  noted  that  some  of  these  figures-of-merit,  such 

as  R,  R'  and  Rgt,  apply  to  more  than  one  and  in  the  case  of  R'  to  as 

many  as  four  different  thermal  environments  and  criteria  of  thermal  stress 

resistance.  It  should  be  noted  that  the  present  review  does  not  include  the 
23  2 

f igure-of-merit  , R’ ' ' = E/St  which  implicitly  is  contained  within  the 
figure-of-merit,  R' ' ' ' . For  this  reason,  the  figures-of-merit  presented 
are  non-redundant.  On  the  other  hand,  a figure-of-merit,  R^  is  introduced, 
equal  to  the  reciprocal  of  the  thermal  diffusivity.  Other  factors  being 
constant,  this  figure-of-merit  maximizes  the  time-of -maximum-stress  (i.e., 
time-of-failure)  of  a brittle  material  subjected  to  transient  heating  or 
cooling.  Material  selection  on  the  basis  of  this  parameter  is  approoriate 
to  structures  or  components  which  have  a specific  function  for  a relatively 
short  duration.  Note  also  that  the  velocity  of  sound  in  a material  has 
been  introduced  as  a material  property  which  can  affect  thermal  stress  re- 


sistance. 


For  the  proper  choice  of  figure-of-merit,  it  is  essential  that 
the  heat  transfer  environment,  the  criterion  of  thermal  stress  resistance 
and  failure  mechanism  be  firmly  established,  before  a choice  can  be  made 
of  the  material  with  optimum  thermal  stress  resistance.  This  choice  of 
figure-of-merit  is  rendered  even  more  complex  by  the  temperature  depen- 
dence of  the  physical  properties,  as  well  as  the  temperature  dependence 
of  heat  transfer  mechanisms.  Due  to  these  temperature  dependencies,  the 
relative  thermal  stress  resistance  of  any  two  materials  may  well  be  inter- 
changed for  different  levels  of  temperature  at  which  the  materials  are  sus- 
ceptiple  to  thermal  stress  failure. 

In  summary,  a review  is  presented  of  all  thermal-stress  figures-of- 
raerit  available  at  the  time  of  preparation  of  this  paper.  It  is  hoped 
that  this  review  will  make  a contribution  to  the  general  understanding 
of  the  nature  of  thermal  stress  failure  of  brittle  materials  and  the 
selection  of  the  material  with  the  best  thermal  stress  resistance. 
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4. 


List  of  Symbols 


A.  Material  Properties 

a - coefficient  of  linear  thermal  expansion 

A - constant  in  V = AKj  exp(-Q/RT) 

a - thermal  diffusivity 

c ■*  speed  of  sound 

e - emissivity 

E - Young's  modulus 

n - viscosity 

G - fracture  surface  energy 

K - thermal  conductivity 

Kt  - critical  stress  intensity  factor 

Ic 

X - wave-length  above  which  dielectric  material  is  opaque 

o 

and  below  which  it  is  transparent 
n - constant  in  V = AKj  exp(-Q/RT) 

v - Poisson's  ratio 

Q - activation  energy  for  sub-critical  crack  growth 

Sc  - compressive  strength 

St  - tensile  strength 

B.  Environmental,  geometric  variables  and  physical  constants 

A - cross-sectional  area  of  column 

b - thickness  of  plate,  column,  radius  of  cylinder,  sphere 

or  spherical  cavity 

B - Biot  number,  B = bh/K 

C,C'  - geometric  constants 


fraction  of  total  energy  from  black-body  radiation 
below  frequency, 

maximum  rate  of  internal  heat  generation 
convective  heat  transfer  coefficient 
cross-sectional  moment  of  inertia  of  column 
crack  half-length 
length  of  column,  rod 

number  of  cracks  per  unit  area  of  volume 
maximum  heat  flux  or  flow 


maximum  rate  of  change  of  heat  flux  or  flow 


Stefan-Boltzmann  constant 


stress 


maximum  temperature  of  black-body  radiation 
critical  temperature  difference 
maximum  temperature  difference 


maximum  rate  of  change  of  AT 


t = L/c 
m 
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